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O ■ This is an introductory review about the on-going search for a signal of 

Lorentz Invariance Violation (LIV) in cosmic rays. We first summarise basic 
aspects of cosmic rays, focusing on rays of ultra high energy (UHECRs). We 
discuss the Greisen-Zatsepin-Kuz'min (GZK) energy cutoff for cosmic protons, 
which is predicted due to photopion production in the Cosmic Microwave 
£T) \ Background (CMB). This is a process of modest energy in the proton rest 

frame. It can be investigated to a high precision in the laboratory, if Lorentz 
transformations apply even at factors 7 ~ O(10 n ). For heavier nuclei the 
energy attenuation is even faster due to photo-disintegration, again if this 
process is Lorentz invariant. Bence the viability of Lorentz symmetry up to 
tremendous 7-factors — far beyond accelerator tests — is a central issue. 

Next we comment on conceptual aspects of Lorentz Invariance and the 
possibility of its spontaneous breaking. This could lead to slightly particle 
^ 1 dependent "Maximal Attainable Velocities" . We discuss their effect in decays, 

Cerenkov radiation, the GZK cutoff and neutrino oscillation in cosmic rays. 

We also review the search for LIV in cosmic 7-rays. For multi TeV 7-rays 
we possibly encounter another puzzle related to the transparency of the CMB, 
similar to the GZK cutoff, due to electron/positron creation and subsequent 
inverse Compton scattering. The photons emitted in a Gamma Ray Burst 
occur at lower energies, but their very long path provides access to informa- 
tion not far from the Planck scale. We discuss conceivable non-linear photon 
dispersions based on non-commutative geometry or effective approaches. 

No LIV has been observed so far. Bowever, even extremely tiny LIV effects 
could change the predictions for cosmic ray physics drastically. 

An Appendix is devoted to the recent hypothesis by the Pierre Auger 
Collaboration, which identifies nearby Active Galactic Nuclei — or objects 
next to them — as probable UBECR sources. 
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1 Present address: Institut fur Theoretische Physik, Universitat Regensburg, 
D-93040 Regensburg, Germany. E-Mail: bietenho@ifh.de 
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Preface : This overview is designed for non-experts who are interested 
in a cursory look at this exciting and lively field of research. We intend 
to sketch phenomenological highlights and theoretical concepts in an en- 
tertaining and self-contained form (as far as possible), avoiding technical 
details, and without claiming to be complete or rigorous. 
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1 Cosmic rays 



1.1 Discovery and basic properties of cosmic rays 

Cosmic rays consist of particles and nuclei of cosmic origin with high 
energy. Part of the literature restricts this term to electrically charged 
objects, but in this review we adapt the broader definition, which also 
embraces cosmic neutrinos and photons. 

The discovery of cosmic rays dates back to the beginning of the 20 th 
century. At that time electroscopes and electrometers were developed to 
a point which enabled the reliable measurement of ionising 7-radiation. 
Around 1909 the natural radioactivity on the surface of the Earth was 
already well-explored, and scientists turned their attention to the radia- 
tion above ground. If its origin was solely radioactivity in the Earth, a 
rapid reduction would have been predicted with increasing hight. 

In 1910 first observations on balloons led to contradictory results: K. 
Bergwitz did report such a rapid decrease [I], whereas A. Gockel could 
not confirm it [2J. Still in 1910 Th. Wulf — who had constructed an 
improved electrometer — performed measurements on the Eiffel tower. 
He found some reduction of the 7-radiation compared to the ground, but 
this effect was clearly weaker than he had expected [3j. Although he 
admitted that the iron masses of the tower could affect the outcome, he 
concluded that either the absorption in the air is weaker than expected, 
or that there are significant sources above ground. 

In 1912 V.F. Hess evaluated his measurements during seven balloon 
journeys pf]: he found only a very mild reduction of the radiation (hardly 
10 %) for heights around 1000 m. Rising further he observed a slight 
increase, so that the radiation intensity around 2000 was very similar 
that that on ground. As he went even higher (up to 5350 m) he found a 
significant increase, which was confirmed in independent balloon-borne 
experiments by W. Kolhorster in 1913 up to 6300 m [5]. 

In Ref. [6] V.F. Hess presented a detailed analysis of his measurements 
in heights of (1000 . . . 2000) m, which were most reliable. These results 
agreed essentially with those by A. Gockel, who had also anticipated the 
observation of an intensified radiation when he rose even higher [2] . Hess 
ruled out the hypothesis of sources mainly inside the Earth. He added 
that sources in the air would require e.g. a RaC density in high sectors 
of the atmosphere, which exceeds the density measured near ground by 
about a factor of ~ 20. He concluded that a large fraction of the "pen- 
etrating radiation" does apparently originate outside the Earth and its 
atmosphere. In a specific test during an eclipse in 1912, and in the com- 
parison of data that he obtained at day and at night, he did not observe 
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relevant changes. From this he inferred further that the sun is unlikely 
to be a significant source of this radiation [3] J! which we now denote as 
cosmic rays (the earlier German term was "Hohenstrahlung" ) . 

Another seminal observation was made by P. Auger in the Alps (1937) 
[7]. Several Geiger counters, which were separated by tens or hundreds 
of meters, detected an event at practically the same time. Such a corre- 
lation, in excess of accidental coincidence, was also reported in the same 
year by a collaboration in Berlin [8]. Later B. Rossi commented that 
he had made similar observations already in 1934. P. Auger gave the 
correct interpretation of this phenomenon, namely the occurrence of ex- 
tended air showers caused by cosmic rays: a high energy primary particle 
of cosmic origin triggers a large cascade in the terrestrial atmosphere [7] . 
Based on the comparison of air shower measurements at sea level and on 
Jungfraujoch (in Switzerland, 3500 m above sea level), and some simple 
assumptions about the shower propagation, P. Auger et al. conjectured 
that primary particles occur at energies of at least 10 15 eV [9]. 

Today we know that the primary particles carry energies of about 
E ?a (10 9 . . . 10 20 ) eV. The emerging cascades typically involve (in their 
maxima) 0(1) secondary particle for each GeV of the primary energy, so 
an air shower can include up to O(10 n ) particles. The energy dependence 
of the flux is illustrated in Figure [D Its intensity falls off with a power 
law oc E~ 3 (new fits shift the exponent to —2.7), up to minor deviations 
("knee"[f], "ankle"). We will come back extensively to the behaviour in 
the highest energy sector. 

The composition of the rays depends on the energy. Here we are 
particularly interested in ultra high energy cosmic rays (UHECRs), with 
energies E > 10 18 ' 5 eV, which are mostly of extragalactic origin. There 
one traditionally assumed about 90 % of the primary particles to be pro- 
tons, about 9% helium nuclei (a), plus small contributions of heavier 
nuclei. At lower energies these fractions change. In that regime, the 
magnetic field of ~ 3 /iG in our galaxy confines charged particles as dif- 
fuse radiation, so that the galactic component dominates [12] • In the 
range (10 15 . . . 10 17 ) eV the rays mainly consist of heavier nuclei [13], and 
for (10 12 . . . 10 15 ) eV there are ~ 50 % protons, ~ 25 % a and ~ 13 % 
C, N and O nuclei, according to Ref. [10]. As we consider still lower 
energies, also leptons and photons contribute significantly, see Sections 
2 and 3. We quote here typical numbers from the literature, but many 
aspects of the energy dependent composition are controversial |12j . 

2 Radiation from the sun occurs at low energy, up to O(10 8 ) eV. It won't be 
addressed here, and we adapt a notion of cosmic rays starting at higher energy. 

3 In addition, a small "second knee" has been observed around 10 17 5 eV. It is not 
visible in Figure [TJ but we will comment on it in Subsection 1.3. 
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Figure 1: The flux intensity of cosmic rays as a function of the energy 
(plot adapted from the HiRes Collaboration Jffl$). Over a broad energy 
interval it falls off essentially like flux oc E~ 3 (or oc E~ 2 - 7 , according 
to recent high precision measurements). We are most interested in the 
UHECRs, with energies E > IO 18,5 eV. Around E = 6 • 10 19 eV an abrupt 
flux reduction is predicted. This is the GZK cutoff, to be discussed in 
Subsection 1.3 and beyond. 
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The arrival directions of charged rays are essentially isotropic. This 
is explained by their deviation in interstellar magnetic fields: such fields 
occur at magnitudes of [iG inside the galaxies, and of nG extragalacticly, 
which is — over a large distance — sufficient for a sizable deflection. As 
a consequence the origin of charged rays from far distances can hardly 
be located^ 

The origin of cosmic rays is still mysterious. The first proposal was 
the Fermi mechanism ("second order"), which is based on particle colli- 
sions in an interstellar magnetic cloud [14J. Collisions which are (almost) 
head-on are statistically favoured and lead to acceleration. A later ver- 
sion of the Fermi mechanism ( "first order" ) refers to shock waves in the 
remnant after a supernova [15J. However, these mechanisms can explain 
the cosmic rays at best in part; they cannot provide sufficient energies 
for UHECRs. In particular the first order mechanism could only attain 
about E = O(10 14 ) eV, and it predicts a flux oc E~ 2 — but the observed 
flux is close to a behaviour oc E~ 3 , as Figure [1] shows. A variety of sce- 
narios has been suggested later, for comprehensive reviews we refer to 
Refs. [TB]. They can roughly be divided into two classes: 

• Bottom-up scenarios : Certain celestial objects are equipped with 
some mechanism to accelerate particles to these tremendous ener- 
gies (which can exceed the energies reached in terrestrial accelera- 
tors by at least 7 orders of magnitude). The question what these 
objects could be is puzzling. Pulsarsj [17] and quasars^] [18] are 
among the candidates that were considered, and recently Active 



sources, see Appendix A. However, a convincing explanation for 
the acceleration mechanism is outstanding (see Ref. [19] for a re- 
cent discussion, and Ref. [20] for a modern theory). 

4 Recent developments could change this picture for the UHECRs, see Appendix 
A. Of course, the approximately straight path length increases with the energy, so 
that UHECR directions may be conclusive if the source is nearby. As an example, a 
proton of E = 10 20 eV has in our galaxy a Larmor radius, rijkpc] « £'[EeV]/-B[/zG], 
which exceeds the galaxy radius w 15 kpc (1 EeV = 10 18 eV, cf. Appendix B). 

D Pulsars are rotating magnetised neutron stars. 

6 A quasar, or quasi-stellar radio source, is an extremely bright centre of a young 
galaxy. 

7 An Active Galactic Nucleus is the environment of a super-massive black hole in 
the centre of a galaxy (its mass is estimated around (10 6 . . . 10 10 )M Q , where M Q ~ 
2 • 10 30 kg is the solar mass). It absorbs large quantities of matter and emits highly 
energetic particles, in particular photons, electrons and positrons (we repeat that 
the acceleration mechanism is not known). Only a small subset of the galaxies have 
active nuclei. Our galaxy does not belong to them; the nearest AGN are located in 
Centaurus A and Virgo; this will be of importance in Appendix A. 




attention as possible UHECR 
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• Top-down scenarios : The UHECRs originate from the decay or an- 
nihilation of super-heavy particles spread over the Universe. These 
particles were generated in the very early Universe, so here one 
does not need to explain where the energy comes from (see e.g. 
Refs. [211 [22]). One has to explain, however, what kind of parti- 
cles this could be: magnetic monopoles have been advocated [25] . 
while another community refers to exotic candidates called "wim- 
pzillas" ]24j - but there are no experimental signals for any of 
these hypothetical objects. Moreover, this approach predicts a sig- 
nificant UHECR flux with 7 and v as primary particles, which is 
disfavoured by recent observations: in particular Ref. [25] estab- 
lished above 10 19 eV an upper bound of 2 % for the photon flux. 
The negative impact on top-down scenarios is further discussed in 
Refs. [26l [27]. Ref. [28] had reported earlier measurements with 
similar conclusions. 



1.2 The Cosmic Microwave Background 

We proceed to another highlight of the historic development: in 1965 
A. A. Penzias and R.W. Wilson discovered (rather accidentally) the Cos- 
mic Microwave Background (CMB) [29], which gave rise to a break- 
through for the evidence in favour of the Big Bang scenario. It is es- 
timated that the CMB emerged after about 380 000 years. At that time 
the cooling of the plasma reached a point where hydrogen atoms were 
formed. Now photons could scatter off these electrically neutral objects, 
so that the Universe became transparent (for reviews, see Refs. [3D]). 
The CMB that we observe today — after a period of 13.7(2) Gyr of fur- 
ther cooling — obeys to a high accuracy Planck's formula for the energy 
dependent photon density in black body radiation, 



duo it 2 exp(u / (hsT)) — 1 ' 

where to is the photon energy (for h — 1) and &b is Boltzmann's constant. 
This shape, plotted in Figure [2} is in excellent agreement with the most 
precise observation, which was performed by the Cosmic Background Ex- 
plorer (COBE) satellite [31]. Any deviation from the black body formula 
( 11. ID over the wave length range A ~ (0.5 ... 5) mm is below 0.005 % of 
the CMB peak@ The temperature is identified as T = 2.725(1) K [32]. 

8 Sizable deviations seem to occur, however, at much larger wave length, A > 1 m, 
due to an additional radio background, the details of which are little known. 
We add that there is also a Cosmic Neutrino Background, which decoupled already 
about 2 s after the Big Bang, and its present temperature is estimated as < 1.9 K. 
We do not consider it here because it has no significant impact on cosmic rays. 
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Figure 2: The Planck distribution ( f 1. 1\) of the black body radiation at 
the CMB temperature of T = 2.725 K. The CMB energy density agrees 
accurately with this Planck distribution — it is by far the most precise 
Planck radiation that has ever been measured. 



This implies that the mean photon energy and wave length are given by 
(u) ~ 6 ■ eV and (A) ~ 1.9 mmj^l The resulting photon density in 
the Universe amounts to J °° du (dn 7 /du) ~ 410.4(5) cm -3 . Interesting 
further aspects of the CMB — which are, however, not directly relevant 
for our discussion — are reviewed for instance in Ref. [53] . 



1.3 The Greisen-Zatsepin-Kuz'min cutoff 

In the subsequent year the knowledge about cosmic rays and about the 
CMB led to an epoch-making theoretical prediction, which was worked 
out independently by K. Greisen at Cornell University [34], and by G.T. 
Zatsepin and V.A. Kuz'min at the Lebedev Institute [35]. They expected 
the flux of cosmic rays to drop abruptly when the energy exceeds the 
"GZK cutoff 

(1.2) 

The reason is that protons above this energy interact with background 
photons to generate pions. The dominant channel for this photopion 
production follows the scheme 

p + 7 ^ A(1232 MeV) ~* P + 7r °, , - ,- (1-3) 

— > n + n , n — ► p + e + v e . 

9 The term "microwave" usually refers to wave lengths in the range of about 
1 mm ... 10 cm, so it does apply to the CMB. 

10 Since this "cutoff" is not sharp, it is a bit arbitrary where exactly to put it. Eq. 
(|1.2[) gives a value which is roughly averaged over the literature, and which we are 
going to embed into phenomenological and theoretical considerations. 



E GZK « 6 • 10 iy eV 



S 



These two channels cover 99.4% of the A(1232) decays [32] • If even 
more energy than the threshold for this transition is available, excited 
proton states (p*(1440), p*(1520) . . . ) and higher A resonances (A(1600), 
A(1620), A(1700) . . . ) can contribute to the photopion production as 
well. In these cases one may end up with p + n or also with p + 2ir; at 
even higher energy the production of three pions is possible too. 

It is obvious that a proton with energy E p > E GZ k looses energy 
under photopion production, until it drops below the threshold for this 
process^ We are now going to take a somewhat more quantitative look 
at transition (11. 3p (in natural units, H — c = 1). 

We denote the proton energy E p at the threshold for the A (1232) 
resonance as E Q . In these considerations one refers to the Friedmann- 
Robert son- Walker metrics (see Appendix B) — which is co-moving with 
the expanding Universe — as the "laboratory frame" (we will call it the 
"FRW laboratory frame" )@ We first consider the relativistic invariant 

s = (E + u;) 2 -(p p + p 7 ) 2 

head— on 

= El - f v + 2E lu - 2p p ■ p~ t ^> m\ + AE uj = m\ . (1.4) 

p p and p 7 are the 3-momenta of the proton and CMB photon. We arrived 
at the term 4E uj by assuming a head-on collision (and E 3> m p ). The 
last term refers to a A baryon at rest, which sets the energy threshold. 
Thus we obtain an expression for E q (lv). As an exceptionally high photon 
energy we insert uj = 5(u) = 3 meV — photons with even higher energy 
are very rare due to the exponential decay of the Planck distribution 
( Oil . This leads tc@ 

2 2 

E = A 4uj P L=5M - 5-3 • 10 19 eV . (1.5) 

Further kinematic transformations yield a simple expression for the 
inelasticity factor K, which represents the relative energy loss of the 

11 In the original studies [3U [35] terms like "A resonance" did not occur, but the 
authors knew effectively about the photopion production. They inferred this spectac- 
ular conclusion, although Ref. [34] only consists of 2 pages without any formula, and 
Ref. [35] is just a little more extensive. Its original Russian version has been trans- 
lated but only few people have read (the page number is often quoted incorrectly). 
Nevertheless both works are of course top-cited on a renowned level. 

12 In terms of a CMB decomposition in spherical harmonics, the monopole contri- 
bution sets the temperature of 2.725(1) K. The dipole term is frame dependent; the 
requirement to make it vanish singles out the maximally isotropic frame [32] . which 
we could also refer to in this context. 

13 In eq. (| 1 . 2[) we increased this energy threshold slightly, which is consistent with 
the unlikelihood of exact head-on collisions. 
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proton under photopion production [55] . 



K:=^L= 1 -h- m '- mi ). (1.6) 
E p 2\ s J K J 

This relative loss still refers to the FRW laboratory frame, but formula 
( 11.61) holds for general proton-photon scattering angles. 

Let us now change the perspective and consider the Mandelstam vari- 
able s in the rest frame of the proton, where we denote the photon 4- 
momentum as (u/,p^), 

s = (m p + to') 2 - p' 2 = m 2 p + 2m p cu' . (1.7) 

The photon energy u' is related to uj by the Doppler effect, 

uj' = •juj (1 — fpcosi?) . (1.8) 

v p and $ are the proton velocity and the scattering angle in the FRW 
laboratory frame. The Lorentz factor 7 can take a remarkable magnitude; 
for instance at the proton threshold energy E p = E it amounts to 

7 = ^\e p =e « 10 11 . (1.9) 

VTLp 

By averaging uj' over the angle $ one arrives at [36] 

E 

cj~180MeV--^. (1.10) 

Hence a proton with E p > Eq perceives the CMB photons as quite ener- 
getic 7-radiation0 

Combining eqs. (11.61) . (11.7j) and (ll.lOp we can determine the inelas- 
ticity factor at a given proton energy (and averaged scattering angle), 

1 / mi — ml 



m = jfi- m r is ) 

2 V m p (m p + 2u) / 

0.15 at E p = E , Q = 180 MeV ( . 

0.20 at E P = 2E , cu = 300 MeV . ^ ' ' 

We see that the inelasticity is significant already at the energy threshold, 
and beyond it (gradually) rises further. 



14 The velocity of the proton in the centre-of-mass frame is given by uj' j (uj' + m p ). 
At E p « E it is small on the relativistic scale, hence the proton rest frame referred 
to above is not that far from the centre-of-mass frame. 
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We are now prepared to tackle the question which ultimately matters 
for the GZK cutoff: if a proton with E p > E travels through the Uni- 
verse, how long is the decay time of its energy ? 

This question was analysed by F.W. Stecker (36], who derived the fol- 
lowing formula for the energy decay time r, 

Tk-^L^^ 1 ^-^ 2 ^ • (ii2) 

Uq is the photon threshold energy for photopion production in the rest 
frame of a proton (which has FRW laboratory energy E p ). The product 
of Q with the logarithm emerges from the Planck distribution (11. ip after 
integration. is the total cross-section for the photopion produc- 

tion}^ which was explored experimentally already in the 1950's, so Refs. 
[SUSHI [36] could refer to it. The corresponding experiment with protons 
at rest, exposed to a 7-ray beam of about 200 MeV, had been reported 
for instance in Ref . [37] . The cross-section varies mildly with the energy, 
in the magnitude of cr ~ 0.1 mb. Its profile was reproduced in Ref. [36J, 
along with the measured inelasticity factor K(uj). It is remarkable that 
this simple experiment at modest energy provides relevant information 
about the fate of UHECRs. Based on the knowledge about a, Figure 2 in 
Ref. [SB] displays the mean free path length for the proton, for instance 

^10£ ^ 10M P C - ( L13 ) 

A later fit to the experimental data, which refers directly to the product 
a(u))K(u)), is worked out in Ref. [38 J . It includes corrections like the 
higher photopion production channels, which still reduce the free path 
length a little@ 

Refs. [2U EH] pointed out already that the corresponding energy at- 
tenuation for heavier nuclei is even stronger. The main reason here is 
the photo- disintegration into lighter nuclei due to interactions with CMB 
photons (see Ref. [39] for a detailed analysis). The fragments carry lower 
energies; the energy per nucleon remains approximately constant. Hence 
we are in fact dealing with an absolute limitation for the energy of cosmic 
rays consisting of nucleons0 Numerous works have later reconsidered 

15 Strictly speaking it is cleaner to set the lower bound in this integral to zero 
and rely on the strong suppression below u>o due to <j(ui), but the form (| 1 . 1 2|) is 
more intuitive. Moreover eq. ()1 . 12() simplifies the energy attenuation to a continuous 
process — its discrete nature gives rise to minor corrections. 

16 The same is true for the somewhat higher CMB temperature in scattering pro- 
cesses long ago (although this effect is marginal). Various CMB temperatures have 
been considered in Ref. [35] . 

17 Regarding non-nuclear rays, we stress that neutrinos are not subject to any the- 
oretical energy limit in the CMB, but no UHECR neutrinos have been observed yet. 
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the attenuation of protons and heavier nuclei in the CMB in detail, see 
e.g. Refs. [101 E] an d references therein. 
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Figure 3: The energy attenuation for protonic UHECRs due to photo- 
pion production with CMB photons (plot adapted from the Pierre Auger 
Collaboration U^). Note that 1 EeV = 10 18 eV (see Appendix B). 

At last one might object that protons may still travel very long dis- 
tances with energies above Eqzk if they start at a much higher energy. 
However, at E p ^> E GZ k the attenuation is strongly intensified. (For 
instance the value of the lower integral bound lj in eq. (11.121) decreases; 
once the integral captures the peak in Figure [21 many more photons can 
contribute to the photopion production). As a result, extremely high 
energies decrease rapidly, and the final super-GZK path does not exceed 
the corresponding path length at E p « E drastically. This feature is 
shown in a plot by the Pierre Auger Collaboration, which we reproduce 
in Figure [31 

It seems natural to assume a homogeneous distribution of UHECR 
sources in the Universe, both for bottom-up and for top-down scenarios 
(cf. Subsection 1.1). Then one could expect a pile- up of primary par- 
ticles just below E G zk- As a correction one should take into account 
another process, which dominates somewhat below Eqzk, and which en- 
ables the CMB to reduce the proton energy further: in the energy range 
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5 • 10 17 eV < E < Eqzk the electron/positron pair production 

p + -f^p + e + + e~ (1.14) 

causes this effect, though the inelasticity is much smaller than in the case 
of photopion production (consider the ratio in eq. (11.111) ). 

Ref. [JT] presents an updated overview, which incorporates all relevant 
effects. This leads to Figure H] for the attenuation lengths for various 
primary nuclei. For the proton it is composed as t~ x ~ t~l + f" 1 , which 
refers to photopion production and / e + pair creation (at even lower 
energy adding also the inverse attenuation length due to the redshift 
becomes significant). 




1 10 100 1000 1 10 100 1000 

E [EeV] E [EeV] 



Figure 4: Illustration of the energy attenuation of UHECRs due to pho- 
topion production with CMB photons for protons (on the left), and due 
to photo- disintegration of heavier nuclei (on the right). In both cases also 
e~ I e + pair production is considered (upper curves at high energy). The 
plot on the left even captures the redshift, which is, however, not really 
important at the ultra high energies that we are interested in. These plots 
are adapted from Ref. [J^l]. 

All in all, one concludes the following: if protons — or heavier nuclei 
- hit our atmosphere with energies E > Egzk, they are expected to 
come from a distance of maximally 



100 Mpc . (1.15) 



If we rise the "cutoff" to 10 20 eV (cf. footnote [10]) the range decreases 
to £ max ~ 50 Mpc. For instance this distance reaches out to the Virgo 
cluster of galaxies (its centre is about 20 Mpc from here). £ max is a 
large distance compared to our galaxy (or Milky Way) — our galactic 
plane has a diameter of about 30 kpc. But £ max is short compared to the 
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(co-moving) radius of the visible Universe. On that scale the sources of 
super-GZK radiation should be nearby, 

(-RMilky Way — 15 kpc) <C £ max <C (-Ryisible Universe ~ 14 Gpc) . (1.16) 

Despite the e + /e~ pair production, some pile-up seems to occur at 
energies E<E GZK (cf. footnote [22]). The status of the subsequent dip is 
discussed in Refs. [13]. As we mentioned in footnote[3J, a "second knee" in 
the cosmic flux was observed by various collaborations around 10 17 ' 5 eV. 
As a possible interpretation it could be the pile-up at the lower end of 
the pair creation threshold [12J. 

However, in this review we want to focus on the issue of super-GZK 
energies in cosmic rays. Not even in our vicinity — given by the range 
(I1.15P — we know of any acceleration mechanism, which could come near 
such tremendous energies. Therefore the detection of super-GZK events 
could represent a mystery — and perhaps a hint for new physics. So let 
us now summarise the actual observations in this respect. 

1.4 Observations of super-GZK events 

Even before the GZK cutoff was established from the theoretical side, one 
spectacular event had been reported by the Volcano Ranch Observatory 
in the desert of New Mexico [2] . The energy of the cosmic particle that 
triggered this event was estimated at 10 20 eV, which exceeds already the 
GZK cutoff (11.21) . Refs. [3U[35] both quoted this report and commented 
that they do not expect any events at even higher energies (Greisen added 
that he found even that event "surprising"). In fact, identifying the 
energy of the primary particle is a delicate issue — which we will address 
below — and the corresponding techniques were at an early stage. 

In any case this initiated a long-standing and controversial challenge 
to verify the validity of the GZK cutoff. In 1971 yet another super-GZK 
event was reported from an air shower detected near Tokyo [45j@ This 
inspired a large-scale installation in Akeno (170 km west of Tokyo), which 
is known as AG AS A (Akeno Giant Air Shower Array); a historic account 
is given in Ref. |46j. 

Meanwhile also the Fly's Eye project in Utah went into operation, 
and in 1991 it reported an event with 3 • 10 20 eV primary particle energy 
[47j . Up to Ref. [IS], this is the highest cosmic ray energy ever reported. 
In macroscopic units it corresponds to 48 J f^l 

18 For that event the energy was estimated even as 4 • 10 21 eV, but it is not often 
quoted in the literature. 

19 It is popular to compare this to the kinetic energy of a tennis balls. Its mass is 
typically 57.6 g, so that these 48 J are the kinetic energy for a speed of 147 km/h, 
which is somewhat below the speed attained in a professional game. 
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Up to the 21 st century numerous new super-GZK events were de- 
tected, in particular by AGASA. According to the data of this collab- 
oration, the spectrum — i.e. the flux of UHECR rays — continues to 
follow the power-like decrease in the energy that we saw in Figure [U 
with hardly any extra suppression as Eqzk is exceeded. This contradic- 
tion to the theoretical expectation triggered an avalanche of speculations. 
The AGASA results have been considered essentially consistent with the 
data of further air shower detectors in Yakutsk (Russia, at Lake Baikal) 
[48] and in Haverah Park (England, near Leeds) [49]. An overview of 
their super-GZK statistics is given in Figure [51 
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Figure 5: The statistics of GZK and super-GZK events reported by the 
observatories AGASA, Haverah Park and Yakutsk, against the energy in 
units of 10 19 eV (Figure adapted from Ref. I5$J ). 

This observation disagrees, however, with the data of the HiRes (High 
Resolution Fly's Eye) observatory, which concludes that the cosmic rays 
do respect the GZK cutoff [STJ. HiRes is the successor experiment of 
Fly's Eye since 1994, with a refined technology!^*! 

A comparison of the cosmic ray spectra as determined by AGASA 
(1990 - 2004) and by HiRes (1997 - 2006) is shown in Figured Com- 
mented overviews over the results of these collaborations are included 
in Ref. [401 E2]- Figure also shows that the overall statistics of these 

20 Initially the HiRes results seemed to agree as well, but here and in the following 
we refer to the results presented later by this Collaboration. 
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Figure 6: A direct comparison of the cosmic flux at ultra high energy as 
measured by the Observatories AG ASA (upper symbols) and HiRes (lower 
symbols). The plot on the left shows the flux as a function of the energy E. 
On the right it is rescaled by a factor E 3 ( cf. Figure Q]], which renders it 
approximately constant up to Eqzk- We see that the continuation beyond 
Egzk is controversial (plots adapted from Ref. JEEj)- 

UHECR observations is modest. Note that the flux above 10 12 eV is 
around 10 primary particles per minute and m 2 , but above 10 18 ' 5 eV 
it drops to 0(1) particle/(km 2 year) (as indicated in Figured]), so the 
search for UHECR takes patience. One may question if there is really a 
significant contradiction between HiRes and the other groups mentioned 
above. In fact, an analysis in Ref. [51] concludes that this discrepancy 
between AGASA and HiRes might also be explained by statistical fluc- 
tuations, if the energy scale is corrected in each case by 15 % (which is 
well below the systematic uncertainty). 

Nevertheless a large community assumed that a contradiction is likely 
and wondered about possible reasons. An obvious suspicion is that this 
may actually be a discrepancy between different methods. 

• To sketch these two methods, we first address the air shower, which 
is illustrated in Figure CD At an early stage the huge energy trans- 
fer of an UHECR primary particle onto the molecules in the atmo- 
sphere generates — among other effects — a large number of light 
mesons (pions, kaons . . . ), which rapidly undergo leptonic (or pho- 
tonic) decays. Of specific interest for the observer are the muons 
emerging in this way: thanks to the strong time dilation they often 
survive the path of about 25 km to the surface of the Earth (al- 
though their mean life time is only 2- 10~ 6 s). Another consequence 
of the extremely high speed is that the air shower is confined to a 
narrow cone — typically the directions of motion of the secondary 
particles deviate by less than 1° from the primary particle direction 
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in the shower maximum (a few hundred meters from the core). The 
method applied by AGASA, Yakutsk and Haverah Park is an ar- 
ray of numerous Cerenkov detectors — spread over a large area on 
ground — which identify the trajectories and energies of highly en- 
ergetic secondary particles. They use tanks with pure water, where 
Cerenkov radiation is amplified by photomultipliers. Detecting in 
particular a set of leptons fi ± , e ± (and photons) belonging to the 
same shower — with precise time of arrival, speed and direction 
- provides valuable information about the shower. Its evolution 
in the atmosphere, and ultimately the primary particle energy, are 
reconstructed by means of sophisticated numerical methods. (Of 
course, the information recorded on Earth is still insufficient for a 
unique reconstruction, so that maximal likelihood methods^] have 
to be applied to trace back the probable scenario.) 

• The method used by HiRes - - and previously by Fly's Eye - 
observes the showers in the atmosphere, before they arrive at Earth. 
The showers excite N2 molecules of the air, which subsequently 
emit UV or bluish light when decaying to their ground state (with 
wave lengths A ~ (220. . .440) nm). Although this light is weak, 
it is visible from Earth to powerful telescopes, along with light 
collecting mirrors, at least in nights with hardly any moonshine or 
clouds. (These telescopes are structured in a form similar to the 
compound eye of an insect.) In 1976 Volcano Ranch achieved the 
first successful observation, which was confirmed by the detection 
of the same air shower on ground. This air fluorescence 
is emitted isotropically. It is also illustrated in Figure El along 
with the Cerenkov radiation in the air as a further effect. The 
direct nature of this observation is a valuable virtue, but an obvious 
disadvantage is the limitation of the detection time to ~ 10 %. In 
this method, the observable height grows with the energy of the 
primary particle, which complicates the interpretation of the data 
[52]. 

We add a couple of qualitative remarks, without going into details. An 
obvious question is how protons as primary particles can be distinguished 
from heavier nuclei. For the latter the air showers arrive at the maximal 
number of particles at a higher point. The depth -D max in the atmosphere 

21 A synopsis of this method is given in Ref. 32]. 

22 We adapt this term from the literature, although it actually suggests that the 
primary particle must be a photon. The appropriate term for this effect would be 
"scintillation" , or more generally "luminescence" . 
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Figure 7: The evolution of an air shower, as illustrated by the HiRes 
Collaboration FlCtf . We recognise the fluorescence light (UV or blue), the 
Cerenkov radiation in the air, and the production of light mesons, which 
rapidly decay into leptons and photons (see zoom on the left). 

to that point can be parametrised as [55j 127] 

[Anax (proton) - D max (A)) oc hi A , (1.17) 

where A is the atomic mass number. This point can approximately be 
identified by either of these two methods: for the fluorescence method 
this observation is quite direct, and on ground it is characterised by the 
electron/muon ratio. Further criteria — such as the time profile of the 
signal and the curvature of the shower front — are reviewed in detail in 
Refs. [56J, and summarised in Ref. [T2j . 

In particular the Fly's Eye record of 3 • 10 20 eV seemed to originate 
from a quite heavy primary nucleus, such as oxygen. However, the iden- 
tification is not easy at all, and in practice the criteria are not always 
consistent. Only cosmic 7-rays can be distinguished quite clearly from 
other primary particles (they penetrate much deeper into the atmosphere, 
hence the air shower maximum is closer to Earth). 

In the course of a long flight through the CMB, heavy nuclei are ex- 
pected to break apart (photo-disintegration), as we mentioned before. 
Therefore a dominance of protons in the UHECRs suggests that the pri- 
mary particles come from relatively far distances, and vice versa. 
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The status of UHECR observations to this point is reviewed in detail 
in Refs. [TBI |4"U] . In light of the dilemma between the results obtained 
with these two techniques, the Pierre Auger Collaboration designed a 
new project in Argentina (near Malagiie, province of Mendoza) with the 
goal to clarify the situation [42] . Its planning started in 1992 and it 
is in stable operation since January 2004, while part of the equipment 
has still been installed. The concept of the Pierre Auger project is the 
combination of both techniques described above: 

• On Earth it involves 1600 detectors with 12 tonnes of water, where 
three photomultipliers monitor the Cerenkov light caused by air 
showers. The installation of these tanks has been terminated in 
2008. They are distributed over an area of 3000 km 2 on a triangular 
grid of spacing 1.5 km (for comparison: air showers due to UHECRs 
have diameters of about 6 km at the terrestrial surface). 

• In addition 24 telescopes are searching for fluorescence light from 
4 well-separated sites — all of them have already been operating 
since 2004. 

The Cerenkov array is suitable for collecting large statistics, while the 
telescopes are important for a reliable calibration of the energy measure- 
ment. This can be achieved best based on the data from air showers, 
which are observed in both ways, hybrid events, by evaluating a vari- 
ety of correlations [57J. With a single detection technique the energy 
calibration has been notoriously problematic. 

In Refs. [571 ESI EH] the Pierre Auger Collaboration presented data 
collected until the middle of 2007, which we reproduce in Figure El The 
exposure to this point already exceeds the total exposure which was ac- 
cumulated by HiRes and AGASA by about a factor of 2 resp. 4. The 
systematic uncertainty in the energy measurement is estimated around 
22 % and the statistical error around 6 % (which is rather harmless in this 
context). The vertical arrows in this plot mark the upper limits for the 
energy bins with 84% C.L. [TT], based on the Feldman-Cousins method 
for Poisson distributions [59]. The small labels indicate the number of 
events detected in the corresponding bin, and the error bars represent 
the statistical uncertainty of the flux. Below E GZK (as given in eq. (11.21) ) 
we recognise the power-like behaviour of the energy spectrum. In the 
intervals 4 • 10 18 . . . 4 • 10 19 eV, and beyond 4 • 10 19 eV its slope was mea- 
sured oc E~ 2 m ^ resp. oc E~ A ' 2 ^ 2 \ which rules out a single power law by 
6 standard deviations. Hence for E > Eqzk the flux drops clearly be- 
low the extrapolated power law. Nevertheless a considerable number of 
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Figure 8: The UHECR spectrum according to the data by the Pierre 
Auger Collaboration collected from January 2004 until August 2007 (plots 
adapted from Ref. J77)/). The lower panel reveals a clear reduction of the 
rescaled flux beyond Eqik, roughly consistent with the HiRes data and 
the predicted GZK cutoff. Nevertheless there is a considerable number of 
new super-GZK events (the upper panel displays the number of events). 
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new super-GZK events have been observed^ Confronted with the fluxes 
measured earlier by AGASA and HiRes, these new data are closer to the 
latter. They are certainly consistent with an extra damping once the 
energy threshold for the A resonance is exceeded. On the other hand, 
the sizable number of super-GZK cosmic rays asks for an explanation 
and keeps the door open for speculations. An interesting interpretation 
of the UHECR arrival directions was published by the Pierre Auger Col- 
laboration in November 2007, see Appendix A. 

In this review we focus on Lorentz Invariance Violation (LIV) as one 
attempt to explain a possible excess of UHECRs compared to the theo- 
retical prediction. We should add, however, that the "boring" outcome of 
full consistency with the GZK cutoff is by no means disproved — recently 
that scenario has actually been boosted, see again Appendix A. 



2 Lorentz Invariance and its possible vio- 
lation 

So far we have assumed Lorentz Invariance (LI) to hold. If fact, it played 
a key role in the derivation of the GZK cutoff: the transformation of the 
scattering process of an UHECR proton and a CMB photon to the rest 
frame of the proton established the link to the experimentally known 
cross-section. Also the determination of the inelasticity factor for the 
proton under photopion production depends on LI, since it refers again 
to the proton rest frame (regardless whether one relies on theory, as in 
eq. (11.111) . or on laboratory experiments). The same holds for the photo- 
disintegration of heavier nuclei, and for electron/positron pair creation. 

LI was therefore crucial for the analysis of the energy attenuation of 
super-GZK cosmic rays. Here we relied on Lorentz transformations with 
extreme Lorentz factors of the magnitude 

7 = Z 10 11 , (2.1) 



v 



see eq. ( 11. 91) (we still set c = 1). There are no direct experimental tests 
if LI still holds without any modifications for such extreme boosts; ac- 
celerator experiments are limited to Lorentz factors 7<O(10 5 ). If the 
validity of the GZK cutoff will ultimately be confirmed, then this ob- 
servation may be considered an indirect piece of evidence for LI under 
tremendous boosts. If, on the other hand, the final analysis reveals a 



23 If one rescales the flux with the factor E 3 — as in Figure[6] — the shape turns into 
a small peak below Eqzk (which could be interpreted as a pile-up, as we mentioned 
in Subsection 1.3), and the suppression beyond the threshold appears weaker. 
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mysterious excess of super-GZK events, it might indicate new physics. 
In that case, one point to question about the standard picture is LI - 
so let us take a closer look at it. 

LI is a central characteristic of relativity: 

• In Special Relativity Theory (SRT) LI holds as a global symmetry. 
The transformation rules were derived by H.A. Lorentz (1899 and 
1904), H. Poincare declared them a Law of Nature, A. Einstein 
provided a consistent physical picture (1905) and H. Minkowski 
embedded it into the geometry of a 4d space (1907). 

• In General Relativity Theory (GRT) LI holds as a local symmetry 
(A. Einstein, around 1915). In each space-time point the frame 
can be chosen so that locally the metrics takes the Minkowski form 
g = diag (1, — 1, —1, —1), which holds globally in SRT. 

It is therefore tempting — and possible — to write GRT in the 
terminology of a gauge theory (see for instance Ref. |60j, Section 
6.4). However, as a quantum field theory this formulation is not 
renormalisable: the renormalisation group flow does not lead to a 
conformal field theory (i.e. to scale invariance) at high energy, as 
reviewed recently in Ref. [61]. Therefore our established description 
of particle physics (local quantum field theory) is incompatible with 
GRT, which applies to very low energy and describes long-range 
gravity successfully. 

Hence LI plays an essential role in both, SRT and GRT, although 
this role is not the same. In this review we mostly address high parti- 
cle energies up to the maximum that has been observed, hence SRT is 
in general the appropriate framework. At specific points, however, we 
are also going to comment on conceivable connections to GRT — this 
becomes relevant in view of speculations up to the Planck scale P^l 

In general it is untypical for global symmetries to hold exactly, so 
in SRT one could puzzle why this should be the case for LI (and the 
related CPT invariance, see Subsection 2.1)@ On the other hand, gauge 
symmetries are exact, hence in this respect the GRT picture appears 
helpful (as long as no violation of LI is observed). 

Nevertheless for our further discussion the framework of SRT is essen- 
tial, since it allows us to apply particle physics as described by quantum 

24 For approaches to quantum gravity, which discuss in particular the role of LI and 
its possible violation, see e.g. Refs. [521 [551 [Ml 155] . 

25 In this spirit, the Standard Model of particle physics appears more natural in 
the now established form, which incorporates neutrino masses, and therefore only 
approximate chiral symmetry. 
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field theory. In field theory, the manifestation of LI as a global symme- 
try is well established. Some field $(x), which may represent a scalar, 
a spinor, a 4-vector or some tensor, transforms in a (finite dimensional) 
representation D of the Lorentz group 5*0(1,3), 



see for instance Ref. [HD] ■ In particular scalar quantities — like the La- 
grangian of a relativistic field theory — are Lorentz invariant (-D(A) = 1). 

In contrast to the rotation group, the Lorentz group is non-compact. 
The extrapolation of LI to arbitrarily high energy leads to the generic 
UV divergences in quantum field theories. The formulation on non- 
commutative spaces — to be addressed in Subsection 3.3 — was orig- 
inally motivated as an attempt to weaken these UV divergences. Further 
early suggestions to deviate from LI in quantum field theories include 
Refs. ESI. 



According to the CPT Theorem, particle physics is invariant if three 

discrete transformations are performed simultaneously: 

C : flips particles into the corresponding anti-particles and vice versa. 

This changes the signs of all couplings to gauge fields, hence C is 

denoted as charge conjugation. 
P : parity, space reflection, x — ► — x . 
T : inversion of the time direction, t — > — t . 

The CPT Theorem states that, assuming locality and LI, CPT invariance 
is inevitable. This was recognised first by W. Pauli and G. Liiders 
and put on a rigorous basis by R. Jost [58] . 

The basic idea is an analytic continuation of the Lorentz group to imag- 
inary time (Wick rotation), so one arrives at a reflection at the origin of 
4d Euclidean space. This reflection captures P and T, and the Wick ro- 
tation implies that C is included as well. The combination amounts to an 
anti-unitary symmetry transformation. Applying the CPT transforma- 
tion twice, one also obtains a derivation of the Spin-Statistics Theorem. 
(A historically and mathematically precise account is given in Ref. [69]). 

Only recently O.W. Greenberg turned the consideration around and 
proved the following Theorem [70J: Any breaking of CPT invariance nec- 
essarily entails a Lorentz Invariance Violation (LIV). 

This Theorem is very general, it applies whenever CPT invariance 
breaks in any Wightman function. 



$(x) 



f/(A) t $(x)f/(A) = J D(A)$(A- 1 a;) 
U unitary , A G 50(1,3) , 



(2.2) 




2.1 Link to the CPT Theorem 
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On the other hand, if we only assume LIV, then CPT symmetry may 
be broken or not, i. e. there are many ways to break LI while leaving CPT 
symmetry intact. The leading LIV terms are linear in the energy resp. 
momenta, and these terms break CPT too; they are CPT odd. If we only 
admit LIV terms starting quadratically in the energy, then CPT can be 
preserved [7T], so it is natural to assume 

LIV|cpt odd OC E , LIV|cPT even OC E 2 . (2.3) 

2.2 Status of experimental LI tests 

According to Greenberg's Theorem, LI tests can be performed indirectly 
by verifying CPT invariance. Experiments probe in particular the CPT 
prediction of identical masses for particles and their anti-particles. This 
holds for instance for e + , e~ up to 8- 10~ 9 , and for p, p up to 1CT 8 (in both 
cases this is the bound for the relative deviation with 90% C.L.) [32J. 
On-going tests investigate the strangeness driven oscillation Kq <-> Kq 
and conclude [73] 

\M Ko -M Ko \< 5.1 ■ 10~ 10 eV (95 % C.L.) , (2.4) 

which corresponds to a remarkable relative precision of O(10~ 18 ). 

Similar high precision LI tests via CPT measure the magnetic mo- 
ments of elementary particles and anti-particles [71]. A further indirect 
method to test LI in part deals with the equivalence principle [75] . How- 
ever, direct tests are even more powerful, and we will explain below why 
we are more interested in CPT- conserving LIV parameters. 

The status of numerous experimental LI tests by means of a variety 
of methods has been reviewed extensively in Refs. [721 [7HJ [77] . Atomic 
physics gives access to excellent precisions for specific LIV parameters. 
Such parameters describe, for instance, a particle spin coupling to a con- 
ceivable "tensor background field". We will come back to these objects in 
the following Subsection. Here we anticipate that a large number of such 
parameters can be introduced as additional couplings in the Standard 
Model. In light of the discussion in Subsection 2.1, we stress that these 
LIV terms are local, but they can be even or odd under CPT transfor- 
mation. According to Ref. [78] they are all bounded by O(10 -27 ), in a 
form which factors out the Planck scale (12. 6p . 

These precisions are impressive, and all results are in agreement with 
LI — no LIV has ever been observed. One may wonder if even higher 
precision is still of interest, or if it is already sufficient for all phenomeno- 
logical purposes. In the following we would like to argue that this is not 
necessarily the case. The precision achieved so far may be insufficient 
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with respect to scenarios, which could drastically affect physics on the 
Planck sale, or even the propagation of cosmic rays. 

• We first address the CPT conserving class of LIVs, which exclude 
terms of 0(E) [71] . so we can assume the violation effect to be 
oc E 2 , as in relation (I2.3p . Let us consider the case where LI is 
violated dramatically on the Planck scale, where one expects both, 
gravity effects and quantum effects to be strong. For instance the 
magnitude of (dimensionless) LIV effects at such energies could be 

LIV(£ = M Planck ) = 0(1) , (2.5) 

where Mpi anc k is the Planck mass, i.e. the energy scale set by the 
gravitational constant G, 

M Planck = 1/y/G ~ 1.2 • 10 28 eV ~ 22 . (2.6) 

Accelerator experiments reach energies E < O(10 13 ) eV, so they 
could only be confronted with such LIV effects < O(10~ 30 ). It is 
therefore conceivable that a CPT conserving LIV has been over- 
looked up to now, although it plays an essential role on the Planck 
scale. 

• We move on to the case where LI and CPT are violated. Then 
we expect LIV effects oc E. Let us assume that even a surpris- 
ingly strong LIV parameter of this kind, say with a magnitude of 
O(10~ 23 ), has been overlooked. It would be amplified to O(10 -8 ) 
on the Planck scale, so even there it would be of minor importance. 

Therefore the CPT conserving LIV terms are more interesting. In ad- 
dition to this estimate of magnitudes one might also argue that breaking 
one fundamental law of standard quantum field theory — which has been 
a principal pillar of physics since the 20 th century — is a very daring step 
already. Hence it appears reasonable to study this step in detail, rather 
than damaging yet another highly established principle on top of it. 

Cosmic rays represent indeed a unique opportunity to observe effects 
not too far below the Planck scale. The point is not only that they carry 
the highest particle energies in the Universe, but — even more impor- 
tantly — their extremely long flight, which could accumulate tiny effects 
of new physics also at low or moderate energy (we will sketch examples in 
Subsection 3.4). If there is any hope to find experimental hints related to 
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theoretical approaches like string theory! 26 ! or (loop) quantum gravity] 27 ! 
then it is likely to involve cosmic rays (although also LHC might have a 
chance |85])o These hypothetical theories tend to install new fields in 
the vacuum, in addition to the Higgs field of the Standard Model. We 
are now going to glimpse at the mechanism how such new background 
fields could lead to spontaneous LIV. 



2.3 Standard Model Extension (SME) 

A systematic approach to add local LIV terms to the Standard Model 
has been worked out by A. Kostelecky and collaborators starting 10 years 
ago [87]. This collaboration provides a kind of encyclopedia of such 
terms, which is reviewed — along with its experimental bounds — in Ref. 
[88] . The original motivation emerged from string theory [89], which is, 
however, not needed for the resulting Standard Model Extension (SME). 

We take a look at a prototype of a LIV term in this extension. Con- 
sider (in a short-hand notation) the Lagrangian of some fermion field 
ip(x), which may represent a quark or a lepton, 

C = i-fa^ - g$W> - ig'G^Yd^ + ... (2.7) 

We start with the free kinetic term, followed by the Yukawa coupling g 
to the Higgs field 0. The Standard Model assumes to undergo Sponta- 
neous Symmetry Breaking (SSB), so that one of its components takes a 
non-vanishing expectation value, say (0o) > 0. This leads to the fermion 
mass g(<po)- A mass term cannot be written into the Lagrangian directly 

26 String theories start from a higher dimension and first need a spontaneous LIV in 
that space to explain why we experience only 4 extended dimensions. This requires a 
non-perturbative mechanism, which is not well understood. One attempt to formulate 
string theory beyond perturbation theory is the lOd IIB matrix model (or IKKT 
model) [79] ; a possible spontaneous splitting into extended and compact directions is 
discussed in Refs. [50], and for its 4d counterpart in Refs. [51] . 

As a different aspect, LIV effects on the Casimir force in a 5d brane world are discussed 
in Ref. [82] . 

27 Conceivable loop quantum gravity effects on cosmic rays are discussed in Refs. 
[83j . So far we can restrict parameters in specific quantum gravity approaches based 
on the bounds on LIV parameters [54"] . 

28 We mention at the side-line that the scattering of UHECR primary particles on 
atmospheric molecules takes place at higher energies in the centre-of-mass frame than 
any scattering in accelerator experiments. In particular the proton-proton collisions 
at LHC (with 7 TeV per beam) will be equivalent in the centre-of-mass energy to the 
scattering of a 10 17 eV proton on a fixed target. This happens about 2.5 • 10 5 times 
per second in the terrestrial atmosphere, in addition to the cosmic rays hitting other 
astronomical bodies, which obviously invalidates popular concerns about the safety 
of LHC [55] . Here, however, we are interested in interactions of UHECRs with CMB 
photons, which involve modest centre-of-mass energies, see Subsection 1.3. 
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if we want to keep the freedom to couple the left- and the right-handed 
chirality components of the fermion field independently to a gauge field; 
this is the situation in the electroweak sector of the Standard Modell^l 
Let us now consider the third term in the Lagrangian (|2.7p . which is 
one of the possible extensions beyond the Standard Model. The Yukawa- 
type coupling g' G IR is another free parameter, which couples the fermion 
to a new background field of the Higgs-type, albeit with a tensor struc- 
ture. Again we assume SSB, which could, for instance, lead to 

(G 00 ) > , (G^) = (otherwise) . (2.8) 

This additional term will obviously distort the fermion dispersion relation 
compared to the Standard Model, which amounts to a spontaneous LIV. 

In contrast to the Higgs field we are now confronted with the ques- 
tion if new background fields — which may have a vector or a tensor 
structure — are Lorentz transformed as well. For a simple change of 
the observer's inertial frame, all fields — including the background fields 
- are transformed, and £ remains invariant. However, the question 
if particles, which are described by fields like if), really perceive LI de- 
pends on a transformation of these fields only in a constant background 
[9T)l 172] . Kostelecky et al. denote this as an "active" (or "particle") 
Lorentz transformation, and in this respect LI may break spontaneously, 




In this way we can construct a non-standard dispersion relation for 
any particle, depending on its coupling to the tensor field G^ y or further 
new background fields. The inclusion of gauge interactions proceeds in 
the familiar way (one promotes some global symmetries to local ones by 
means of covariant derivatives). 

Of course the LIV parameter g' is just one example - - the mass 
term and kinetic term in the Lagrangian (12. 7p could both be extended 
by including a background field for each element of the Clifford algebra 
(further terms are written down in Section 3.5). Another example, which 
was considered earlier [91 J, is the addition of an extra gauge term of the 
Chern-Simons type to modify qedH 

£cs = -\e^p,A u F pa = \B-A, Gv/ = // > 0) . (2.9) 

29 0n the other hand, a fermion mass term —m't/jijj is allowed in a vector theory like 
QCD, where both chirality components of a quark couple to the gluons in the same 
manner. 

30 Here one re-introduces a distinction, which one had hoped to overcome since the 
historic work by Einstein. 

31 In non-Abelian gauge theory, LIV through Chern-Simons-type terms is studied 
in Refs. [92]. 
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Here is a vector of dimension mass (and B = V x A). Each component 
has vanishing covariant derivatives in all frames. This vector introduces a 
"preferred direction" , where also all the ordinary derivatives $„p M vanish 
- Ref. (91] associates it with the preferred frame in a galaxy. It identifies 
an astrophysical bound of \i < O(10~ 33 ) eV based on radio galaxies (the 
method will be sketched in Subsection 3.5). Recently there have been 
attempts to achieve even higher precision for this bound based on the 
rotation of CMB photon polarisation vectors [93J. 

Note that £qs breaks also CPT invariance spontaneously. Examples 
for CPT conserving LIV terms in the fermionic and gauge part of the ex- 
tended Lagrangian are the LIV term in eq. (12. 7p . or K ilvpa F >iV F p ' 7 (where 
K is another background tensor field). Kostelecky et al. have identified 
more than 100 LIV parameters in this way, including CPT breaking terms 
[94]. The resulting model preserves a number of properties of the Stan- 
dard Model, like energy and momentum conservation, gauge invariance 
and locality. For massive fermions at low energy, and small LIV parame- 
ters, also energy positivity and causality are safe [95J. In the framework 
of SRT the photon is identified with the Nambu-Goldstone boson of the 
SSB of Lorentz symmetry. That collaboration also discussed various sce- 
narios for an interpretation in the GRT framework [871 EU [88] . 

In this parametrisation, the current status is summarised in Ref. [78] 
as follows: all (dimensionless) LIV coefficients are < 10~ 27 (as we men- 
tioned before). For those which break CPT the upper bound is as tiny as 
10 -46 , which further motivates our focus on CPT even terms. A detailed 
overview of the various bounds is given in Ref. [77] . 

An obvious question is why LIV should become manifest only at huge 
energies, although SSB is typically a low energy phenomenon. In princi- 
ple even a LIV detected a low energy would not necessarily imply LIV 
at high energies. Of course we are interested in the opposite situation, 
where LIV is significant only at high energies, so we have to assume the 
SSB to persist over all energies that we consider. This setting requires 
tiny LIV parameters multiplying momenta of the fields (as in the exam- 
ples above). Then the deviation from LI is only visible at huge momenta 
and a contradiction to known phenomenology can be avoided. 

Although consistent, this scenario implies a severe fine tuning prob- 
lem. On tree level one might hope for a somehow natural suppression of 
the LIV parameters by a factor m ew /M Planck ~ 10~ 17 (where m ew is the 
electroweak mass scale) [HE] • If one includes 1-loop radiative corrections, 
however, this effect is lost and the magnitude of LIV effects is just given 
by 0(g^ w ) (g cw being an electroweak Yukawa coupling) [97] , which is far 
from the suppression needed. If we consider LIV nevertheless, we have to 
add this problem to the list of unsolved hierarchy problems, like the small 
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vacuum angle 9 in QCD, the small particle masses (on Planck scale) or 
the small cosmo logical constant. On the other hand, the Standard Mod- 
els of particle physics and cosmology are well established, despite the 
presence of these unsolved problems. 

2.4 Maximal Attainable Velocities (MAVs) 

Within the same framework we now switch to the pragmatic perspective 
of Ref . [71] . We saw in the previous Subsection that LIV can be arranged 
for. Now we assume that this mechanism has been at work, and we 
proceed to an effective Lagrangian, which includes explicit LIV terms. 
Such an effective Lagrangian may also capture hadrons as composite 
particles (like Chiral Perturbation Theory). 

Following Ref. [71] we leave, however, other Standard Model proper- 
ties intact as far as possible. We therefore require CPT and gauge in- 
variance to persist, and 50(3) rotation symmetry to hold in a "preferred 
frame" @ Moreover we do not break power counting renormalisability, 
i. e. we only add terms of mass dimension < 4. 

Let us outline which terms are admitted by these conditions. For a 
bosonic field <3> we can add a purely spatial kinetic term 

1 3 

£ efr = • • • + 2 ^ W a K 6 W 6 ) , (e ab = e ba ) . (2.10) 

i=l a,b 

(Equivalently we could add a term with only time derivatives instead). 
For a fermion field we insert a spatial kinetic term as well, where the LIV 
parameters can be distinct for positive or negative chirality, 

£ eff = • • • + kH • d Ml + 75) + £-(1 - 7 5 ) ] V • (2-11) 

Gauge interactions require covariant derivatives also in these non-standard 
terms. For the pure gauge part we first consider the Abelian gauge group 
U(l), where the electric and the magnetic field are given by E % = — E 0t , 
B l = ^e % i k Fji e . The above conditions allow for three independent terms] 33 ! 

E 2 -B 2 , E ■ B , B 2 . (2.12) 

32 In contrast to Ref. [SI], the following considerations on cosmic rays can hardly 
employ the preferred frame in a galaxy, because the primary particle paths extend 
over much larger distances. One might refer to the frame which renders the CMB 
maximally isotropic, cf. footnote \T2\ 

33 Without insisting on CPT invariance, the Chern-Simons term (|2.9|) would be 
added to this list. 
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The first two among these terms are LI, so we work with the third one. 
Also in Yang-Mills gauge theories Ref. [7T] uses the corresponding term 



C cS = --- + g"J2B a B a , (2.13) 

a 

where the index a runs over the generators of the Abelian factors in the 
gauge group. 

In this way one obtains a quasi-Standard Model with 46 LIV param- 
eters. In view of the above list the large number may come as a surprise 
at first sight, but it can be understood by the numerous new fermion gen- 
eration mixings (for one case, the phenomeno logical impact on neutrino 
oscillation will be discussed in Subsection 2.5.4). Coleman and Glashow 
verified that all these parameters preserve the vanishing gauge anomaly, 
so that the model remains gauge invariant on quantum level. 

We proceed to another prototype illustration and consider a neutral 
scalar field [71 J. In some background the renormalised propagator G is 
written as 

- itr 1 = (p 2 - m 2 )f(p 2 ) + efg(p 2 ) . (2.14) 

Geometrically we assume the usual Minkowski space with c = 1, p 2 = 
E 2 — p 2 , and mo is the renormalised mass at e — 0. The functions / and 
g are not specified, but / obeys the conventional normalisation, and g 
is adapted to it, /(mg) = g{rrify = 1 (both functions are smooth in this 
point). 

Let us now turn on a tiny LIV parameter e and consider its leading 
order. The poles in the propagator are shifted to 

771 2 — *2 i 2 ->2 ->2 2 , 2 4 

E ~ p + m — ep ~ p c P + m c P 
with m = , c P = 1 - e (to 0(e)) . (2.15) 

The renormalised mass is modified, itlq — > m. However, we are more 
interested in the feature of the dispersion relation (12.151) : the parameter 
cp takes the role, which is usually assigned to c. From the group velocity 



dE \p\ 



d\p\ 



m 2 c1 



c P (2.16) 



it is obvious that this is the speed that the particle approaches asymp- 
totically at large \p\. 

Following this pattern, each particle type P can pick up its own Max- 
imal Attainable Velocity (MAV) cp; it might slightly deviate from the 
speed c, which establishes the Minkowski metrics. A tiny value of the 
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dimensionless parameter e justifies the above linear approximations, and 
it is compatible with the desired setting where LIVs only become notice- 
able at huge momenta. We mentioned before that this imposes a new 
hierarchy problem (which we have to live with) and that we are referring 
here to an effective approach, which captures composite particles, unlike 
Subsection 2.3. 



2.5 Applications of distinct MAVs 

Let us now address some applications that we are led to if we assume cer- 
tain particles to posses MAVs distinct from c. We still follow Ref. [7lj and 
focus on applications of potential interest in cosmic ray propagation]^! 



2.5.1 Decay at ultra high energy 

We consider the possibility that some particle with index decays into a 
set of particles with indices a = 1, 2 ... . At ultra high energy all masses 
are negligible, but we assume individual MAVs for the particles involved, 
c , c a . We further define c min as the minimal MAV among the decay 
products. Then the energetic decay condition reads 

\P (0) I CO = £ \P W | C a > C min £ \P (a) I > I 
a a 

=^ Co > c m ; n := ™ in C a . (2-17) 

The interesting observation is that tiny differences in the MAVs can be 
arranged such that a usual decay cannot take place anymore if the pri- 
mary particle carries ultra high energy. This happens when the condition 
( I2.17P is violated, i.e. if c is smaller than any of the c a . 

On the other hand, new decays could be allowed at high energy, such 
as the photon decay 7 — * e + + e~ , see Subsections 2.5.2 and 3.5, or the 
radiative muon decay — * e~ + 7. A further example is the inverse 
/3-decay p n + e + + u e , which tightens the bounds on specific LIV 
parameters [99] based on the new results for UHECR that we review in 
Appendix A. 

2.5.2 Vacuum Cerenkov radiation 

Next we consider a particle with a MAV that slightly exceeds the speed 
of light, i.e. it exceeds the MAV of the photon in vacuum, 

^ = l+e>l. (2.18) 
34 It has also been suggested to test MAVs in air showers [98] . 
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This particle can be accelerated to a speed v > c T Then it will emit 
vacuum Cerenkov radiation, so it slows down intensively. The energy 
required for vacuum Cerenkov radiation has to fulfil 

„ m m , 

E > -. ~ —= , 2.19 

y/l " (c 7 /c P ) 2 727' 
where m is the particle mass, and we consider the leading order in e. 

As we discussed in Section 1, we know from cosmic rays that protons 
support an energy up to O(10 20 ) eV over a considerable path, hence this 
energy does not fulfil inequality (12.191) . Coleman and Glashow infer 

2 

771 

e p <^\ E =io^eV-5-10- 23 , (2.20) 

where m p and e p are the mass and the e-parameter (see eq. (12.181) ) of the 
proton. 

Cosmic electrons and positrons are observed only up to E w 1 TeV, 
hence the Cerenkov constraint on their e-parameter is much less strin- 
gent, e e < 10~ 13 . On the other hand, the absence of the photon decay 
7 — > e~ + e + constrains e e from below [100J. In total one arrives at 

_ lo- 5 < ( £i - l) < 1(T 13 . (2.21) 



2.5.3 Impact on the GZK cutoff 

In this generalised framework we now reconsider the head-on collision 
p + 7 -> A (1232) of relation (TO) . For the MAVs we could check all 
kinds of scenarios — they are all conceivable in the effective ansatz of 
Subsection 2.4. Here we pick out the particularly interesting option that 
only the proton has a slightly slower MAV, 

c 7 = c A = l, Cp = 1 - e < 1 , (2.22) 

(in the notation of eq. (12.201) : e = —e p ). The condition for a A resonance 
in the scattering on a CMB photon was written down before in eq. (I1.4p . 
The following requirement still holds, 

m 2 A < (E p + uj) 2 - (p - uj) 2 ~E 2 p -p 2 + 2uj(E p + p) , (2.23) 

where E p , p = \p\ are the proton energy and momentum, and u is the 
photon energy, all in the FRW laboratory frame. The evaluation of 
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the right-hand-side is now modified due to the non-standard dispersion 
relation of the proton, 

El- P 2 {l-ef = ml{l-eY -> E 2 p - p 2 ~ m\ - 2eE 2 p , (2.24) 

where we assumed E p m p and |e| <C 1. If we additionally make use of 
£ p > u, we arrive at a new energy threshold condition, where one term 
is added compared to eq. (11.41) . 

m l - m l + 2e#p < 4:uE p . (2.25) 

Once more, the extra term is suppressed by the factor e, but it can 
become powerful at very large energies. The consequence of this modi- 
fication is illustrated in Figure [9j the right-hand-side increases linearly 
with the proton energy E p and reaches the standard threshold (for e — 0) 
at the value Eq given in eq. ( 11. 51) . A small e > leads to an increased 
threshold energy, and above a critical value e c the A resonance is avoided, 



to ^ 2lo 2 

2 2 l' J= 5 (v)=3 meV 



, , — i . i_ > u,' •> me \ — 3 " 10 . (2.26) 

2E m A — m p 

Hence the critical parameter e c is tiny indeed, and this is all it takes 
to remove the GZK cutoff, at least in view of the dominant channel of 
photopion production^ (This is of course compatible with the condition 
( I2.20p . which only requires e > —5 • 10~ 23 .) The experimental bounds 
for the LIV parameters in the fundamental SME Lagrangian are below 
10~ 27 (cf. Subsection 2.3). However, it is not obvious how this translates 
into bounds on effective MAV differences (remember that here we started 
from a renormalised propag ator)@ 

Let us stress once more that this is not in contradiction to the fact 
that this A resonance is observed at low energy — in that case the masses 
contribute significantly to the energies and the above approximations are 
not valid any longer. 

For e > e c the dominant candidate for a photopion production channel 
is an excited proton state, 

p + 7 -> p*(1435) -> p + n° . (2.27) 

But at ultra high energy p* can only decay if c p > c n , according to 
condition (j2.170 i 37 l With a suitable choice of the MAVs we could close 
that channel too, and so on. 



35 The idea that LIV could invalidate the GZK cutoff was first expressed in Refs. 
[101] , Later it has also been discussed in the SME framework [102J. 

36 A possible transition from SME parameters to hadrons by means of a parton 
model was discussed in Ref. [103] . 

37 Bounds on based on UHECR are discussed in Refs. 65, 104 . 
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Figure 9: An illustration of the threshold for the term AuE p to enable pho- 
topion production through the A(1232) resonance: in the standard situa- 
tion (all MAVs are 1) this threshold is constant, m\ — m^ ~ (807 MeV) 2 . 
For e > the threshold squared picks up a contribution oc E'i, see eq. 
A2.25\) , which requires a higher proton energy, and e > e c ( illustrated by 
the dashed-dotted line) closes this channel for photopion production. 

Just before Ref. [7T] appeared, Ref. [TB] attracted attention to the ob- 
servation that the top 5 super-GZK events known at that time all arrived 
from the direction of a quasar (cf. footnote [6]). This inspired Coleman 
and Glashow to suggest an unconventional hypothesis: the super-GZK 
primary particles could be neutrons. Assuming particle specific MAVs, 
the following arguments yield a consistent picture (see also Ref. [105J ) : 

• For c n < c m i n (where c m in refers to the /3-decay products in the 
sense of relation (12.171) ) the decay of the neutron can be avoided. 

• Since we already assumed a slightly reduced neutron MAV, we 
might add c n < ca to protect the neutron from the GZK cutoff. 

• Unlike the proton, the neutron is hardly deflect by interstellar mag- 
netic fields, hence the scenario of an approximately straight path 

- which is required to make sense of the observation in Ref. [18] 
for quasars at large distances — becomes realistic. 

This suggestion is currently not discussed anymore, but its consis- 
tency is beautiful enough to be reviewed nevertheless. Nowadays 0(100) 
super-GZK events have been reported (see Section 1), and the quasar 
hypothesis is out of fashion as well. However, the hypothesis of clustered 
UHECR directions has recently been revitalised, see Appendix A. Also 
the idea of electrically neutral primary particles is still considered, see 
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e.g. Refs. |1U6] I^I They include pions 7r° [53], photons |1U7] (though Ref. 
[25] puts a narrow bound on their UHECR flux), neutrinos [1U8] and 
hadrons of (u, d, s) quarks, which could be stable in SUSY theories }109] . 

2.5.4 Impact on neutrino oscillation 

If the MAVs can deviate from c, there are three bases for the neutrino 
states: 

• the flavour basis with the states \u e ), \u T ), 

• the basis of the neutrino masses, 

• the basis of the MAVs. 

In this framework neutrino oscillation could occur even at vanishing 
neutrino masses, due to a flavour mixing in the MAV basisj^l However, 
in light of experimental data this scenario has soon been discarded [11 1L 
112J, so we start here from the usual point of view that the observed 
neutrino oscillation is evidence for non-vanishing neutrino masses. 

We do consider, however, the possibility that the mass driven oscil- 
lation could still be modified as a sub-leading effect by additional MAV 
mixing. This scenario was investigated |113j based on the data by the 
MACRO Collaboration in Gran Sasso (Italy). The study can be simpli- 
fied by focusing on the oscillation <-> i/ T , which is most relevant for 
the observation of cosmic neutrinos. In this 2-flavour picture the two 
parameters 

Av : = |MAV(z/i) -MAV(i/ 2 )| 
9 V := mixing angle of and \v T ) in the MAV basis (2.28) 

modify the life time of v^. 

As an example, Figure [10] shows the survival probability of over 
a distance of 10 7 m (roughly a path across the Earth passing near the 
centre) at Av = 2 ■ 10~ 25 (still in natural units, c = 1). The three curves 
refer to the MAV mixing angles with sin26 l „ = 0, ±1. We see that 
the effect of such a mixing is very significant for neutrino energies of 
0(100) GeV. If we assume neutrino masses m< 1 eV, this corresponds 
to a Lorentz factor 7 > 10 11 , similar to the primary proton in a UHECR 
(see eq. fll.9p ). Cosmic neutrinos in this energy range occur @ and their 
behaviour is again suitable for a sensitive test of LI at ultra high 7-factors. 

38 However, neutral primary particles are not required for the hypothesis discussed 
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Av = 2 10^, sin2 6 v = ±T| 




Figure 10: The survival probability for a neutrino over a distance of 
10 7 m, for a MAY difference Av = 2 ■ 1(T 25 (see eq. MTM) ). The plot 
is adapted from Ref. 1113\/ . It shows the curves for mixing angles 9 V = 
(only mixing due to the neutrino masses) and 9 V = ±7r/4 (red dashed 
line and green dotted line), where 9 V parameterises a mixing with v T in 
the MAY basis. We see a strong effect on neutrinos of O(10 n ) eV. 

The MACRO Collaboration detected upward directed muons gener- 
ated by neutrinos in the state through the scattering on a nucleon, 

+ n fi~ + p or + p — > fi + + n . (2.29) 

By monitoring multi-Coulomb scattering they reconstructed the energy 
E^ and finally E v . They observed 58 events with E v > 130 GeV and 
compared this flux to the one of at low energies. It turned out that the 
oscillation curve obtained solely by mass mixing (corresponding to the 
curve for sin 29 v = in Figure [TOl) works very well. Including the addi- 
tional parameters (12.281) does not improve the fits to the data. Therefore 
the scenario of additional «-> v T mixing due to different MAVs is not 
supported, and Ref. [113] concludes 

Av < 6 • 10~ 24 . (2.30) 

Of course the exact bound on Av depends on the mixing angle — it 
becomes even more stringent when | sin2#„| approaches 1. 



in Appendix A, since it assumes nearby UHECR sources. 

39 A discussion in the SME terminology is given in Ref. |110j . 

40 One suspects that there are cosmic neutrinos at much higher energies as well — 
they are not restricted by any GZK-type energy cutoff (cf. footnote [TTJ) nor deviated 
by magnetic fields, so they could open new prospects to explore the Universe. Several 
projects are going to search for them systematically, see Section 4. 
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This result is in agreement with similar analyses based on the Super- 
Kamiokande K2K data [112[ 1114] . Ref . |112] introduced a parameter n 
to distinguish scenarios for the origin of neutrino mixing as follows: 



The fit to the K2K data led to n = —0.9(4), supporting the standard 
picture and strongly constraining the alternatives. Similarly Ref. |114j 
generalised the scenario corresponding to n = 1 by admitting vector and 
tensor background fields — as in Subsection 2.3 — and arrived qualita- 
tively at the same conclusion. 

Regarding the SME mixing parameters, a phenomenological test with 
v e V/i data has been presented in Ref. |115j . On the theoretical side 
the Vn «-> v T oscillation of Majorana neutrinos is discussed in Ref. |116j . 

3 Cosmic 7-rays 

In this section we are going to consider the photons themselves as they 
appear in cosmic rays — so far they entered our discussion only in the 
background. Their identification as primary particles is on relatively safe 
grounds (cf. Section 1). 

3.1 Another puzzle for high energy cosmic rays ? 

The highest energies in cosmic 7-rays were observed around _E 7 > 50 TeV 
(see e.g. Refs. |117j ). They originate from the Crab Nebula, which is the 
remnant of a supernova — recorded by Chinese and Arab astronomers 
in the year 1054 — at a distance of 2 kpc. 

The sources of the strongest 7-rays reaching us from outside our 
galaxy are blazars. Blazars subset of the AGN, see footnote mR a 
few hundred blazars are known. A prominent example is named Markar- 
ian 501; in 1997 the High Energy Gamma Ray Astronomy (HEGRA) 
satellite detected 7-rays with E 1 > 20 TeV from an outburst that this 
blazar had emitted |118j . It is located at a distance of 157 Mpc — this 
can be determined from the redshift (z ~ 0.034). Here also the direction 
of the source is known, in contrast to far-travelling charged cosmic rays 
(gravitational deflections are small). Another example, which is notewor- 
thy in this context, is Markarian 421, at a distance of 110 Mpc, which 
emitted photons that arrived with ~ 10 TeV |119j . 

41 Appendix A addresses the hypothesis that AGN could emit UHECRs as well. 




-1 mass mixing (standard) 

energy independent mixing 

1 LIV resp. violated equivalence principle. 



(2.31) 
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It was suggested that the observation of these highly energetic 7-rays 
could pose a new puzzle, which is similar to the GZK cutoff [120J. Again 
the question is why the CMB can be transparent for rays of such high 
energies. In this case, one expects electron/positron pair creation of an 
UV photon in the ray and an IR photon in the background, 

7uv (£) + 7irH -> e + + e- , (3.1) 

if sufficient energy is involved. This pair creation is followed by inverse 
Compton scattering on further CMB photons, triggering a cascade. In 
the centre-of-mass frame both photons in the transition (13.1 j) have the 
energy 

E . . 

uj = — = ory , (3.2) 

7 

where 7 is the Lorentz factor between centre-of-mass and FRW labora- 
tory frame (cf. Section 1). The condition for pair creation is obvious, 

uj 2 = Eu > m 2 e . (3.3) 

If we insert once more the exceptionally large CMB photon energy that 
we referred to before in Sections 1 and 2, uj 0.003 eV ~ 5 (u) , the 
threshold for E is almost 10 TeV. If we repeat the considerations about 
the CMB photon density and the cross-section for this pair creation, 
in analogy to Subsection 1.3, also the observed energies E of photons 
coming from far distances (outside our galaxy) might be puzzling [121] . 
The mean free path length £ 77 is given by a formula similar to eq. f 1 1 . 1 2 j) , 

, dnjdu [ 4Eu> , . . , oA . 

aw 2 — / dssa^{s) . (3.4) 



Eyy(E) 8E 2 Juj min (E) U J(2m e ) 

As in eq. (11.41) . s is one of the (LI) Mandelstam variables. The differential 
photon density of the CMB follows Planck's formula (II. ip . The minimal 
CMB photon energy is u min (E) = m\j 'E , see eq. (13. 3p . The boundaries 
for the angular integral over s are also obvious, and a(s) is the cross- 
section. It peaks just above the threshold s = (2m e ) 2 , and at high 
energy (s ^> ml) it decays as [21] 

1 s 

a ^ K 7 ln 2^2 ' ( 3 - 5 ) 



s 



Ref. |122j evaluated this formula and found for instance for E = 10 TeV 
a path length £ 77 w 10 Mpc, so that the radiation from Markarian 421 
and 501 appears amazing. Ref. [122] suggested a generalisation of for- 
mula (13. 4p to a LIV form, inserting ad hoc a deformed photon dispersion 
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relation, which extends £ JT However, the question if this puzzle persists 
in the LI form, after taking all corrections into account, is controversial 
|123[ I124[ [76] . One issue is our poor knowledge about the radio back- 
ground radiation, which dominates over the CMB at wave lengths >lm 
(cf. footnote [8]). Ref. [21] studied the effective penetration length Cascade 
of the photon cascade under different assumptions for the extragalactic 
magnetic field and the radio background. Over a broad interval for these 
two parameters, Figure 14 in that work implies 

100 Mpc at E = 30 TeV , . 

300 Mpc at E = 10 TeV . ^ ' 

In light of these numbers the observations of Markarian 421 and 501 
appear less puzzling, but the issue is still under investigation. As in the 
case of the GZK cutoff we discuss possible solutions for such a puzzle - 
if it exists — based on LIV. 

Also this multi- TeV 7 -puzzle could in principle be solved by assuming 
different MAVs (cf. Sections 2.4, 2.5) for the two particle types involved 

c e = c 7 + e e . (3.7) 
For a head-on collision the condition (13.31) is modified to 

Euj > m 2 e + e e E 2 + 0(e 2 e ) , (3.8) 

in analogy to eq. (12.251) . Therefore a parameter e e > could increase the 
energy threshold and — if e e reaches a critical value e c — exclude the 
electron-positron pair creation. Then the Universe becomes transparent 
for photons of higher energy. At the energies that we considered above, 
this happens for 

£e > £c = 2 — -|_E=20TeV, ^=5M=0.003eV — 2-10 15 . (3.9) 

This critical value is much larger than the one for the proton in eq. (I2.26p . 
However, what matters is that e c is well below the upper bound on e e 
that one obtains from the absence of vacuum Cerenkov radiation, given 
in eq. (12.211) . Therefore the MAV scenario could provide a conceivable 
solution to the (possible) multi-TEV 7-puzzle as well. 

Ref. |126j contains a discussion in the spirit of SME and finds that 
LIV, in addition to modifying the minimal energy threshold, could also 
imply an upper energy bound for this pair production. 
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3.2 7-Ray- Bursts (GRBs) 

We now consider cosmic 7-rays of lower energy and focus on j-Ray- 
Bursts (GRBs). These bursts are emitted in powerful eruptions over 
short periods, usually of a few seconds or minutes (in extreme cases 
of (0.1 . . . 10 3 ) s). They typically include photons with energies _E 7 « 
(10 4 . . . 10 8 ) eV. Temporarily they are the brightest spots in the sky. 
The size of these sources is rather small. About the GRB origin we 
only have speculations, such as the merger of two neutron stars, or of 
two black holes, or of a neutron star and a black hole, or hypernovae (a 
theoretically predicted type of supernova, where an exceptionally mas- 
sive star collapses). There are, however, sophisticated models for the 
GRB generation, dealing in particular with fireball shock scenarios or a 
gravitational core collapse 

GRBs have been observed since 1967. Most of them were discovered 
first by satellites (such as Hubble); subsequently ground-based telescopes 
could focus on them and identify the distance to the source. Today also 
the initial detection from Earth is possible, in particular from the obser- 
vatory in La Silla (Chile) |130j . which enables an immediate evaluation of 
the redshift. A spectacular GRB was observed in 2005 from the satellite 
SWIFT: it came from a distance of ~ 4 Gpc; these photons travelled to 
us directly from the Universe when it was only 9 • 10 8 years old. 

The very long journey of photons with different energies but similar 
emission times is suggestive for a test of the photon dispersion relation. 
This test probes if the group velocity t> 7 is really independent of the 
momentum |120j . The simplest approach for a deviation of this standard 
assumption is the introduction of a photon mass m 7 > (which requires 
a modification of the Higgs mechanism). This obviously leads to 

P 

Vy{p) = ^ const. (3.10) 

y/p z + 

as sketched in Figure [TT] (c = 1 is still the photon MAV). From the obser- 
vation that the photons in a GRB arrive at approximately the same time 
one can deduce m 7 < 10~ 6 eV |120j . However, the bound obtained in 
the laboratories is superior by many orders of magnitude, m 7 < 10~ 17 eV 
[32]. 

Therefore the option of introducing a photon mass is not attractive in 
this context — even if it exists below the experimental bound, its effect 
would be invisible to us in GRBs. In the following we stay with m 7 = 0. 
The MAV ansatz E 2 = p 2 c 2 does not help as a theoretical basis for a pos- 
sible momentum dependence of the speed of photons. On the effective 



42 It has also been proposed that the sources of GRB and UHECR could coincide 
P2SEE9]. 
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Figure 11: Qualitative comparison of various dispersion relations in the 
IR regime: standard photon (mass zero in commutative space), massive 
photon (in commutative space), massless scalar particle in an NC space 
and massless photon in an NC space. The latter corresponds to the 1-loop 
result, but its feature changes non-perturbatively, see Subsection 3.3.2. 

level one can still study LIV modified Maxwell equations and evaluate the 
dispersion together with the rotation of the polarisation plane [1311 1132] , 
see Subsection 3.5. However, we now proceed to another theoretical ap- 
proach for new physics, including a non-linear photon dispersion relation, 
which is highly fashionable. That approach is basically different from the 
effective action concepts that we considered in Section 2@ 

3.3 Non- commutative (NC) field theory 

As a brief introduction, we embed this booming field of research into a 
historic perspective. Galilei's system of coordinate transformation did 
not have any constant of Nature (hence it was in some sense more rel- 
ativistic than the relativity theory of the 20 th century). Einstein did 
introduce the speed of light c as such an invariant; it is a corner stone 
of both, Special and General Relativity Theory. Now there are specu- 
lations that there might be yet another constant of Nature in addition 
to c. Ref. |134j denotes this line of thought as "Doubly Special Relativ- 
ity" (according to Ref. |133j it emerged as an effective model for "varying 
speed of light" in a quantum space-time). This second constant is usually 
assumed to be a large, observer-independent energy. 

One often refers to gravity effects becoming important in field theory 

43 A further alternative are "varying speed of light theories" , which are reviewed in 
Ref. [Ti3]. 
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at very large energy. This is supposed to be a generic property which 
is likely to hold even though we do not have an established theory of 
quantum gravity [61]. Therefore it appears natural to relate this second 
constant to the Planck mass (12. 6p (though there may be alternatives 
[135] ). Geometrically this means the absence of sharp points, but a space 
where points are washed out on the scale of the Planck length 

^pianck = Mp^ nck ~ 1.6 • 1(T 35 m . (3.11) 

(J. von Neumann dubbed it "pointless geometry", though this term ap- 
plies to the phase space in standard quantum mechanics as well). 

It can be compared to the lattice regularisations, which truncates 
short distances ad hoc. For numerous quantum field theories this is the 
only regular isat ion scheme, which works beyond perturbation theoryl 44 ! 
it is not only a basis for the numerical measurement of observables, but it 
actually defines these theories on the non-perturbative level. Nevertheless 
hardly anyone would assume a lattice structure to be physical — it could 
at best mimic the geometry on the Planck scale. A more subtle way to 
introduce a geometric UV cutoff places field theory on a space with a non- 
commutative geometry. Let us briefly review the historic construction in 
the first publication on this subject |138] . 

We know that quantum mechanics is endowed with a natural discreti- 
sation of angular momenta based on the uncertainty relation. The idea 
is now to make use of this mechanism for space coordinates as well: the 
coordinate components should be described by operators with discrete 
spectra, as it is familiar for the angular momentum components. This 
can be achieved from a higher dimensional perspective. We start from a 
5d Minkowski-type space with the invariant quantity 

S — Xq ^3 ^4 * ^3.12^ 

We specify a surface S = const. > 0, which is a 4d de Sitter space inside 
the 5d light cone. Now we list the generators of transformations that 
leave S invariant] 45 ! 



h 

L 3 = —(xid2 — £2<9i) , invariant : x\ + x\ , Xq , 23 , 24 , 

Mi = ih(x di — Xido) , invariant : Xq — x\ , x 2 , x 3 , x 4 , 

X\ = -(£i<9 4 — x 4 di) , invariant : x\ + x\ , x Q , x 2 , x 3 , 

T = ia(xod4 — x^dg) , invariant : x% — x\ , x\ , 22 , ^3 • 

44 The "fuzzy sphere" |136j has recently been considered as a possible alternative, 
but numerical work in that regularisation |137j is at an early stage. 

In general we use natural units, but at this point it is instructive to make an 
exception and write h explicitly. 
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The operators Li, L 2} M 2 , M 3 and X 2 , X 3 are analogous. They all con- 
serve 5d LI. Li and Mj generate the familiar rotations and Lorentz boosts, 
so they preserve in addition also 4d LI (at x± — 0). The latter is broken, 
however, by the Xi and T. Here a parameter a enters in a role similar to 
^ in the 4d transformations. 

To show that the spectra of Xi are discrete, we follow the procedure 
which is well-known from quantum mechanics. We write the rotated 
plane in polar coordinates, (xi,x&) = r(sin<y9, cos<^), so that X\ = —iad^. 
For an eigenfunction ip, X±ip = Xip, we make the ansatz 

ip oc exp(iy?A/a) , and ^(v 9 ) = VKv 9 + 2vr) =>- A = na (n E 2Z) . 

As another property that is familiar from the angular momentum algebra, 
these position operators are non- commutative (NC) 

[X h X 3 ] = —e ijk L k , [T, Xi] = —Mi . (3.13) 

That property obviously implies new uncertainty relations for these space- 
time coordinate operators (without involving momenta !) 

AXi AT , AXi AXj > 0(a 2 ) , (i ^ j) . (3.14) 

In this sense a takes indeed the role of a minimal length, and thus of a new 
constant of Nature. A field theory living on this space is necessarily non- 
local over the characteristic extent a. LI is violated in the 4d subspace, 
though its 5d version is preserved. The latter property shows that this 
way to implement a minimal length is far more subtle than a crude space- 
time latticed: some facet of relativity persists, and the mechanism is 
organically embedded into the general framework of quantum physicsF^l 



46 Alternatively, a discrete 4d space-time with Poisson distributed random sites can 
be compatible with local LI [139] . The LI dispersion relation can be reproduced 
exactly, even on regular lattices, if one implements a "perfect lattice action" based on 
the renormalisation group (see e.g. Ref. [140 and refs. therein). 

47 Based on such considerations, it is tempting to speculate that extra dimensions 
could in fact have a physical meaning. The literature on this issue is tremendous; we 
do not try to review it here, but we only mention one point: 

Further arguments for extra dimensions deal with the mass hierarchy problem, i.e. the 
question why the known particle masses are so far below the Planck mass (|2 .6[) (or in 
QCD: m u ^ <C Aqcd)- For instance, fermions naturally undergo a renormalisation, 
which puts their mass non-perturbatively on the cutoff scale (unless some artificial 
fine-tuning of the bare mass is performed) . An extra dimension helps to some extent 
to implement approximate chirality and therefore a small fermion mass [141] . If one 
is willing to assume SUSY in addition, then this mechanism is transmitted to the 
boson masses as well. 
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Snyder constructed the 4-momentum as [138J 



P, = -^ (/i = 0...3), (3.15) 

a £4 

which is consistent with = eyfcpQPj — XjPj), Mi = XjP + TPi. The 
commutators 



[X l ,P j } = ih[6 lJ + (^J (3.16) 



yield Heisenberg uncertainties, which are modified to 0(a 2 PiPj/h). For 
a tiny parameter a we have — on tree level — again the situation that 
the modification is only visible at huge momenta. 

Based on the uncertainty relations (13.141) . the identification of a as 
the Planck length now leads to a consistent picture of the space-time 
uncertainty range as the event horizon of a mini black hole [142] . In 
fact, a (hypothetical) measurement of a length of 0(Lpi anc k) requires 
(according to the Heisenberg uncertainty) an enormous energy density, 
which gives rise to such an event horizon — the notion of detectable 
events then requires inequalities (I3.14p to hold with a ~ £pi anc k- This 
Gedankenexperiment suggests that points should indeed be washed out 
over a range of 0(Lpi anc k). If several directions are involved, as in relation 
(I3.14p . this is practically equivalent to non-commutativity. 

As a variant we may just combine the uncertainty over the gravita- 
tional radius Ax/2 ~ Lp lanck ' with Heisenberg's uncertainty to arrive 
at a modification of the latter, even in one dimension |143j . 

Ax > max(V(2Ap), 2L 2 Planck Ap/h) > L Planck . (3.17) 



However, this modification can lead to dangerous processes, setting in 
at the 1-loop level, where it is not obvious how to avoid a contradiction 
to known phenomenology [144] . We now return to the NC space — al- 
though it is plagued by dangers of this kind as well — keeping in mind 
that basic progress in physics has been repeatedly related to the realisa- 
tion that specific quantities are "unspeakable" |145j . 

Above we followed the historic path to non-commutativity. The cur- 
rently popular version (for a review, see e.g. Ref . 1146] ) inserts for the 
coordinate commutators a constant "tensor field" Jfj 

[X^X V \ = i9 M „ , min(AA^ AX U ) = ^|9J . (3.18) 



48 It is a tensor under a twisted Poincare group, and the action can be called LI if 
one also adapts a twisted particle statistics [147] . 
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Thus we have the case of an "active LIV" (see Subsection 2.3), and agree- 
ment with the general scope of DSR. The question of its CPT invariance 
is discussed in Refs. [HEl UM [T50] . 

The original motivation for introducing NC coordinates was the hope 
to tame the notorious UV divergences in quantum field theory. In fact, 
they are removed in the non-planar diagrams, but in the planar diagrams 
they persist |151j . Moreover, the non-planar diagrams give rise to new 
divergences in the IR sector |152j . which causes the dangerous processes 
that we announced. 

A famous example is the one-particle irreducible 2-point function 
(4>{—p)4>{p)) in the 4d NC A0 4 model. The planar 1-loop contributions 
diverges oc A 2 as in the commutative case (A being a UV cutoff), whereas 
the non-planar diagram picks up an additional phase factor, 

(0(-p)0(p)>l- lo o P , non-planar = ~J J^y fc2 + ~ 2 • (3.19) 

The exponential phase factor is generic for each crossing between lines 
of external propagators and loops. At finite and p v the integral 
J dk^ converges now thanks to the rapid oscillation at large \k p \. As a 
consequence, for this non-planar diagram the UV cutoff turns into |152j 

Als = i/a 2 -?W^ ' (3 ' 20) 

The UV limit A — > oo can in fact be taken in this term, as long as the 
field momentum p remains finite. In addition to the IR divergence due 
to the latter condition, we also see that the limit ||@|| — > is not smooth. 

So far no systematic solution has been found for the treatment of 
multi-loop diagrams, which involve such singularities at both ends of 
the momentum scale. Hence renormalisation does not become easier 
but even more complicated due to non-commutativity. This problem is 
known as UV/IR mixing. Intuitively this is not too surprising in view of 
the uncertainty relation (13.181) . 

Part of the literature deals with a truncated expansion in small ||@||, 
which allows for instance for a generalised formulation of the Standard 
Model [153J. This generalisation may be of interest by itself, but due to 
the UV/IR mixing effects — and the related discontinuous transition to 
standard commutativity, ||0|| ^ — it is fundamentally different from 
complete quantum field theory on a NC space; for instance, locality is 
restored in that truncated expansion. Ref. (154] also refers to such a 
truncated expansion to embed NC field theory into the SME (cf. Sub- 
section 2.3). However, since this alters basic properties, the SME — or 



45 



other low energy effective field theory approaches to include LIV terms 
- and (full-fledged) NC field theory are basically different. 

Despite — or because of — the amazing UV/IR mixing effects, NC 
field theory has attracted a lot of attention; over the last 10 years, about 
2500 works have been written about it. This boom was triggered by the 
observation that it can be mapped onto certain string theories at low 
energy [155J. In that framework one could also relate the non-locality 
range to the string extent; this was the inspiration of Refs. [143] . 

A general analysis of the spatial uncertainty in string theory is given 
in Ref. |156j . and Ref. |157] presents numerical results related to the 
emerging geometry in the IIB matrix model (cf. footnote I2"E|) . 

NC field theory can be viewed as the second main theoretical frame- 
work for LIV, in addition to the approach described in Subsections 2.3 
and 2.4. In both cases it was string theory which inspired a new concept 
in particle physicsF^I but the way that the formulations are elaborated 
now they can be introduced and studied as modifications of standard 
field theory without any reference to strings. 

3.3.1 The photon in a NC world 

Due to non-commutativity, gauge transformations and translations are 
intertwined. As a consequence, even the pure U(l) gauge field picks up a 
Yang-Mills type self-interaction term, i.e. the U(l) gauge field becomes 
NC as well. 

This can be seen easily if we switch to the alternative formulation, 
which uses ordinary coordinates, but all field multiplications are per- 
formed with the star product (or Groenewold-Moyal product) |159] 



Its equivalence to the formulation referred to above can be shown by a 
plane wave decomposition, which allows for a transition between fields 
and Weyl operators (see e.g. Ref. |160j ). The star-commutator [x^, x u }+ = 
iO^ makes this link plausible. In case of a pure U (N) gauge action, 



the commutator term does not vanish, even if the gauge group is U(l), 
since it is promoted to a star-commutator, L4 M , A u ]+ := A^* A v — A u * A^. 

This formulation applies to all NC U (N) gauge theories, but for the 
gauge groups SU (N) it fails: in commutative spaces the algebra of the 

49 The history is reviewed from this perspective in Ref. [158j . 



4>(x) * ij)(x) := 0(x) exp (~ \e^%^(x) 



(3.21) 
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generators A a closes thanks to the identity Tr[A a , A b ] = 0, which does 
not hold, however, for the star-commutator. In view of phenomenology, 
the group U(l) is therefore of primary interest. 

The Yang-Mills type self-interaction in NC spaces leads to a distorted 
photon dispersion relation. On tree level the non-commutativity only 
affects very short distances and therefore the UV behaviour of dispersion 
relations. But on quantum level it causes in addition IR singularities, 
as we pointed out before. Ref. [161] performed a 1-loop calculation of 
the effective potential in NC U(l) gauge theory and derived a dispersion 
relation of the form 

E 2 = p 2 + C^ + 0(g*), 
pop 

pop:= -W Vv ®° Pc > , C = const. (3.22) 

This suggests that a small non-commutativity has a stronger impact than 
a large one, and based on GRB data one could therefore be tempted to 
claim a lower bound |162j F^l which is large compared to the Planck length 
scale, 

||0|| > 10- 44 m 2 ~ (10 13 • L Planck ) 2 • (3.23) 

However, one should better not insist on this apparently spectacular 
inequality (as Ref. [162J also noticed) for a number of reasons: 

• First we stress that this does not rule out the standard commutative 
world, = 0, where no extra term of this kind occurs. Eq. (13.221) 
shows again that the transition ||0|| — > is not smooth. 

Part of the literature deals with an expansion in 1/||0||, which does 
have a smooth limit to a commutative field theory at ||0|| — > oo 

- at least in perturbation theory around the free field vacuum - 
because the extra phase factor exp(ifc At /i!/ p i/ ) (for the crossing of 
propagator lines, see eq. (13.191) ) cancels the non-planar diagrams 
at finite external momentum p. However, that commutative limit 
does not coincide with the standard setting = 0. Moreover, in 
case of SSB the above perturbative argument collapses, because it 
does not refer to an expansion around a vacuum state anymore, 
and therefore NC effects may persist [164|. 

• Such a strong non-commutativity is well above the upper bound 
established from other observations |lb'2j . Does this mean that NC 
field theory is ruled out already ? 

Further aspects make this issue more complicated. 



50 Ref. [163) presents a similar consideration based on 1-loop effects in the potential 
for NC scalars, which becomes strong for very small ||0|| > 0. 
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• A careful evaluation [165] revealed that the NC term in eq. (13.221) 
is actually negative, 

C=~. (3.24) 

Therefore the system appears IR unstable, as illustrated in Figure 
[TTl In string theory there are attempts to interpret such systems 
(a key word is "tachyon condensation" [166j ). but in simple terms 
of quantum physics a system without a ground state does not have 
an obvious meaning. Hence this observation looks like a disaster 
for NC field theory — it seems to be in straight contradiction to 
everyday's experience. 

The behaviour with a positive IR divergence can be realised for 
instance for NC scalar particles. This feature was indeed measured 
numerically in the 3d A0 4 model with a NC spatial plane [167J, 
as sketched in Figure [TTJ A condensation of the finite mode of 
minimal energy then gives rise to a "stripe pattern" formation, i. e. 
some non-uniform periodic order dominates the low energy states, 
which corresponds to SSB of Poincare invariance. In perturbation 
theory the negative IR energy for the NC photon emerges from the 
requirement to fix the gauge. 

There is extensive literature about electrodynamics with star prod- 
ucts on tree level, ignoring this problem. However, the question is 
in which regime non-commutativity effects could be manifest while 
quantum corrections are not visible. Here the surprises on quantum 
level are even more striking than in the case of LIV parameters en- 
tering effective field theory [97] because of the UV/IR mixing. Most 
people who were concerned with this severe problem moved on to 
a supersymmetric version, where the negative IR divergence is can- 
celled by a positive one due to the fermionic partner of the gauge 
field, with an opposite loop correction |165[ 1168] . Does this mean 
that non-commutativity is only sensible if one puts SUSY as yet 
another hypothetical concept on top of it ? 

Even then the corresponding M = 4 super- Yang-Mills theory suf- 
fers from a hierarchy problem for the SUSY breaking scale [149]. 

• From our point of view, the apparent quest for multiple hypotheses 
(NC geometry plus SUSY) to cancel the IR instability makes NC 
field theory less attractive. However, if the 1-loop result has such 
amazing features one may wonder if it should be trusted to extract 
qualitative properties of the system. We therefore reconsidered this 
question on the non-perturbative level [169] . 
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3.3.2 The NC photon revisited non-perturbatively 

A truly non-perturbative investigation of the NC photon — at finite 
non-commutativity and finite self-coupling — has to be carried out nu- 
merically This requires a Euclidean metrics. If we introduce a NC 
Euclidean time, reflexion positivity (as one of the Osterwalder-Schrader 
axioms [1 70 j ) is lost. Then it is not known how to achieve a transition 
to Minkowski signature, so we keep the Euclidean time commutativeF*] 
Moreover we need the (anti-symmetric) tensor O to be invertible on 
the NC subspacej^l hence we only introduce a NC spatial plane (xi,X2) 
with a single non-commutativity parameter 9 

[xi,Xj] = iOeij = 1,2) , (3.25) 

while the plane (£3,0:4) is kept commutative. (Once the Euclidean time 
£4 commutates, this form can always be obtained by means of rotations.) 
This is the only conceptually clean non-perturbative NC setting available, 
and it can be viewed as a minimally NC formulation. 

Next we need a regularisation to a finite number of degrees of freedom 
in order to enable a numerical treatment. As in commutative field theory 
the lattice is suitable for this purpose. Here we review the NC lattice 
formulation briefly, details are explained in Ref. [160j . Although we do 
not have sharp points for the lattice sites at hand, a lattice structure can 
be imposed by requiring the operator identity 

exp (i^-^i) = 4 • (3.26) 

The length a is now the (analog of a) lattice spacing, i.e. the spacing in 
the spectra of the operators X{ (which cannot be diagonalised simultane- 
ously). Again in the spirit of minimal non-commutativity we let the mo- 
mentum components commute (as in Snyder's formulation (13.151) ). Their 
usual periodicity over the Brillouin zone has remarkable consequences, 

\ — e i(ki+2Tr/a)xi e -ikjXj 

= exp (i(ki + —)xi - ikjXj + —)k j [x i , xA 

= ft exp ( — 9(k 2 - kij) => —h E % . (3.27) 
V a J 2a 



51 The opposite case with [xi,t\ = -Xi, [xi,Xj] = |171j is known as n-Minkowski 
space. Work on its impact on dispersion relations is summarised in Ref. [124j . 

52 If this is the case, two operators O at distinct points can be multiplied as |160j 
0(x)0(y) =ir- d (dete)- 1 f d d zO{z) exp[-2i(6- 1 ) A1 „(a;- z)^(y- z)„], where d is the 
dimension of the NC subspace, which has to be even. 
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In the exponent on the right-hand-side we had to keep track of the 
Baker-Campbell-Hausdorff term, which we evaluated by inserting rela- 
tion (13.251) . For fixed parameters 9 and a we infer from this exponential 
factor that the momentum components in the NC lattice are discrete, so 
that this lattice is automatically periodic. Of course, this property is in 
contrast to lattices in a commutative space. Here we inevitably obtain 
an UV and IR regularisation at the same time. 

If our lattice is periodic over N sites in each of the NC directions, 
the momentum components are integer multiples of 2n/(aN), which - 
along with relation (I3.27P — implies 

6 = -Na 2 . (3.28) 

7T 

This shows that the extrapolation of interest is the Double Scaling Limit 
(DSL), which leads simultaneously to the continuum and to infinite vol- 
ume at fixed 9, 

a - (continuum limit) \ at ^ = ^ ^ 



Na — > oo (infinite volume) 

The requirement to take this simultaneous UV and IR limit is again spe- 
cific to NC field theory and related to the UV/IR mixing. Any deviation 
from this balanced condition can lead to a limit, which is not really NC 
(it may lead to 9 — > or oo). This can cause confusion, and in other 
regularisation schemes it is less obvious how to control this equilibrium 
between UV and IR effects. 

In view of lattice simulations, the construction of link variables in the 
NC plane appears difficult — note that a compact link variable should 
now be star-unitary, U * U* = U* * U = 1. However, there is an exact 
mapping [172j onto a U(N) matrix model in one point, the so-called 
twisted Eguchi-Kawai model |173j . which is numerically tractable. 

The DSL of U (1) lattice gauge theory on a NC plane has been studied 
in Refs. |174j . Wilson loops of small areas are real and they follow the 
Gross- Witten area law of commutative 2d U(N — > oo) gauge theory [175] . 
Wilson loops of large area A are shape dependent. For most shapes they 
have a complex phase Aj 9, which corresponds to an Aharonov-Bohm ef- 
fect, if one formally interprets 1/9 as a magnetic field across the plane. 
The same interpretation was suggested long ago based on simple tree 
level considerations |176j . and it is also inherent in the famous Seiberg- 
Witten map [155] onto open bosonic string theory. 

Let us return to the 4d U(l) gauge theory and its formulation on a 
regular lattice in the commutative plane, and a lattice as constructed 
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above in the NC plane. In each direction there is (approximately) the 
same number N of lattice sites. 

In the NC plane the self-interaction and the gauge transformations 
are again non-local. As a consequence, there are gauge invariant open 
Wilson lines, in contrast to gauge theory in commutative spaces. Due to 
the separation of their end-points they carry non-vanishing momenta, 
so these open Wilson lines are suitable order parameters to probe if 
translation and rotation symmetry is spontaneously broken (we refer to 
the discrete translation and rotation symmetry on the lattice, resp. the 
continuous version in the DSL). In Mean Field Approximation, active 
Poincare symmetry will always appear broken because of the constant 
tensor 6^. However, the question whether SSB really sets in or not is 
dynamical. We verified it by measuring the expectation values of open 
Wilson lines. 

Numerical simulations at varying gauge coupling and different lattice 
sizes N revealed that a number of observables — Wilson loops and lines 
and their correlators, matrix spectra as extent of a dynamically generated 
space — stabilise (to a good approximation) for [169J 



This defines the scale for the lattice spacing, up to a proportionality 
constant, which corresponds to the choice of 9. Now the DSL is well- 
defined and it can be explored by increasing N and /3, while keeping 
the ratio N/f3 2 = const. Thus we extrapolate to the continuum and to 
infinite volume at fixed 9. 

Our simulation results led to the phase diagram shown in Figure [12) 
at strong coupling, (3 < 0.35 we find for any iV a symmetric phase, and 
at weak coupling again, but in between (at moderate coupling) there 
is a phase of spontaneously broken Poincare symmetryF^l Apparently 
the strong fluctuations at small (3 avoid SSB on one side, and at weak 
coupling the correlation length exceeds the scale of the mode of minimal 
energy, so that SSB only occurs at moderate coupling. The transition line 
between the weak and moderate coupling roughly follows the behaviour 
(3 C oc N 2 , which means that the DSL curve (I3.30p always ends up in the 
broken phase. So this is the phase that actually describes the NC con- 
tinuum limit. Since it stabilises all observables that could be measured, 
the NC photon could be IR stable after all. 

Figure [13] shows the numerically measured dispersion relations of the 
photon in a NC space. In the symmetric phase at weak coupling it is 

53 This observation agrees qualitatively with later simulations of the 4d twisted 
Eguchi-Kawai model, which were performed with different motivations [177j . 
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Figure 12: The phase diagram for NC QED (pure U(l) gauge theory) 
based on numerical results 11691J . At strong coupling (/3< 0.35), and also 
at weak coupling, translation and rotation symmetry is intact. In between, 
at moderate coupling, there is a phase where it is spontaneously broken. 
Perturbation theory perceives the weak coupling phase. The DSL A3. 29) . 
however, leads to the broken phase, which is therefore physically relevant. 



fully consistent with the IR instability found in perturbation theory, cf. 
Figure HH Hence we agree with that calculation, but it does not capture 
the phase which is really relevant for the NC photon. The limit to a 
continuous, infinite NC space (the DSL) takes it to the broken phase, 
and there Figure [13] shows that the linear dispersion relation is restored. 

This linear dispersion relation is consistent with IR stability and with 
the interpretation of the photon as a Nambu-Goldstone boson of the 
SSB of Poincare symmetry (as it was also suggested in the SME [87]). 
This dispersion relation E(p) was measured in the commutative plane, 
Pi — P2 — 0, p = ps, where the exponential decay of the correlation 
function of open Wilson lines is clearly visible. 

Once IR stability holds, the next step would be to elaborate also 
explicit non-perturbative results for a NC deformation of the dispersion 
in the DSL. Such a study should address finite momentum components 
of the NC directions in the broken phase dispersion relation. This could 
then be confronted with phenomenological data to set bounds on 9. A NC 
deformation in the IR sector could also be highly relevant for cosmic rays 
in view of the background radiation: the GZK cutoff and the multi-TeV 
7-puzzle both depend on the CMB (and radio background) photons|3 



54 Moreover the CMB anisotropy (cf. footnote fT2l and Ref. [33]) might be affected 
by NC effects, if they were sizable in the early Universe |178j . 
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Figure 13: The dispersion relation of the NC photon in the symmetric 
phase (on the left) and in the phase of spontaneously broken Poincare 
symmetry (on the right). The behaviour on the left matches the 1-loop 
result, which is IR unstable, and which we also sketched in FigurelTl\ The 
behaviour on the right was measured in the broken phase, which describes 
an IR stable NC continuum and infinite volume limit {169}. 

The results obtained already [169] suggest the important feature that 
the photon may survive in a NC world after all, without the necessity of 
SUSY to cancel a negative IR divergence. Hence NC QED is an option 
for new physics. Of course, the version with a constant is probably 
still just a mimic of a realistic fuzzy quantum space. 

3.4 Analysis of GRB and blazar flare data 

From the preceding discussion we extract the message that a theoret- 
ical basis for a deformed photon dispersion is conceivable. Similar to 
Subsection 2.4 we now turn to a pragmatic point of view and just as- 
sume some effective deformation. We adapt an obvious parameterisation 
ansatz [120] for the photon energy E at some momentum p, 

M is a very large mass, which emerges somehow, for instance from some 
kind of "quantum gravity foam"@ We keep track of O(EfM), and t> 7 is 
the photon speed (which turns into c = 1 for M — > oo). 

Ref. [181] performed a systematic analysis of the GRB data collected 

55 A general discussion of possible power series corrections to dispersion relations 
and their constraints based on requirements like positivity, causality and rotation 
symmetry is given in Ref. [179] . Ref. [180] studies the impact on the thermodynamics 
of a photon gas. 
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by three satellites named SWIFT, BATSE and HETFj 56 ! and confronted 
the data with ansatz (13.311) . In fact there is a distinction in the time of 
arrival: photons with higher energy tend to arrive earlier (which would 
naively suggest M < 0). However, the analysis has to verify if this could 
also be explained without new physics. We do not know if the photons 
of different energies in a GRB were emitted exactly at the same time (we 
repeat that we have to refer to hypotheses and models for the origin of 
GRBs |127j ). We consider photons in energy bins which differ by AE. 
Without new physics, the delay in the time of arrival can be described 
as 

At (^source (1 

z) (z : redshift) , (3.32) 

where d source is the parameter for the possible time-lag at the source. 
This delay is further enhanced on the way to Earth by the factor (1 + z). 

Then the question is if the data can be fitted better when we add a 
LIV parameter. This leads to the extended parameterisation ansatz 

^source + a Ll\K(z) , d^V : = , J u i (3.33) 



1+z SUU1LB — v--v~, ' — MH 

where H is the Hubble parameter and K(z) is a correction factor, which 
is derived in Ref. [18 lj. Figure [TH (on the left) shows the left-hand- 
side of eq. (I3.33P on the vertical axis, so the question is if the dataEll 
are compatible with a constant, and therefore with LI. Although the 
best linear fit — shown in the plot — has a slight slope, it is obvious 
that a constant behaviour is not excluded. The outcome of a precise 
evaluation is reproduced in Figure dH on the right. The inner (outer) 
depicted area captures the parameters rf SO urce and oliv with 68 % (95 %) 
C.L. The crucial point is that oliv = (absence of LIV) is included in 
the permitted region, and oliv is strongly constrained, which leads to 
the cautious conclusion 

|M| > 1.4 • 10 25 eV « 0.0012 M Planck (3.34) 

with 95 % C.L. In some sense this is a negative result: once again the 
data do not give evidence for new physics. However, it is impressive that 
these observations allow for robust phenomenological insight about an 
energy regime, which is not that far below the Planck scale. 

We add that there are further analyses of GRB data along these lines, 
which do not find evidence for LIV either. In some studies the bound 



56 For instance HETE distinguished for the GRB photons 4 energy bins and detected 
them in time intervals of 6.4 • 10~ 3 s. 

57 The error bars were enhanced by the authors of Ref. |181j compared to the satellite 
data, until they became compatible with a consistent picture. 
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Figure 14: Illustrations of the GRB data analysis in Ref. 11 8 If : on the 
left the left-hand- side of eq. Ii3. 33\) is plotted against the function K(z). 
The question is if the fit requires a non-zero LIV parameter a^iy, i.e. if 
there is a slope with respect to K(z). On the right the window for the 
parameters d SOUTCe (for a lag at the source) and auv is shown for 68% 
resp. 95 % C.L. A strong deviation from the LI scenario auv = is ruled 
out, which implies the bound \3. 3$ (plots adapted from Ref. \181f ). 

for \M\ is set even higher, e.g. Refs. |182[ 1183] analysed the data of sin- 
gle GRBs and concluded \M\ > 0.015 M Planck resp. \M\ > 0.0066 M Planck . 

Very similar investigation are also carried out based on the 7-rays 
emitted in blazar flares. Recently the MAGIC Collaboration used the 
Imaging Atmospheric Cerenkov Telescope to analyse a flare of Markarian 
501 (cf. Subsection 3.1) and arrived even at \M\ > 0.022 M Planck [ 153] . 
The Gamma-ray Large Area Space Telescope (GLAST) [185j is now going 
to search systematically for photons from GRBs and blazar flares in the 
energy interval of 8 keV to 300 GeV, which is promising in view of further 
insight even above these bounds [186J. GLAST has arrived in orbit in 
June 2008. 

If we assume some quantum space-time background, then it is natural 
that it is subject to fluctuations, which may affect the speed of photons 
also stochastically [18TJ . A consequence of this additional effect could be 
the broadening of a flare width, which increases with the distance. 

At last we add that GRBs are also accompanied by a burst of neu- 
trinos, which may lead to additional prospects for LIV tests [188J. 
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3.5 Effective Field Theory and vacuum birefringence 

After discussing NC field theory as a LIV approach in Subsection 3.3, we 
returned to the Effective Field Theory (EFT) framework in Subsection 
3.4. This is the theoretical framework that we had already employed in 
Section 2. It includes the LIV terms, which are not excluded by symme- 
tries, in an energetic hierarchy. This concept works successfully in many 
branches of physics — for instance in Chiral Perturbation Theory — but 
the example of NC field theory shows that it is not without alternative 
(in the latter case, the IR divergences are incompatible with EFT). 

A systematic EFT analysis has to distinguish the mass dimensions 
n of the LIV operators. We summarise this approach briefly, following 
the excellent discussion in Ref. [7S]. We require rotation symmetry and 
gauge invariance, and we introduce a time-like unit vector u. Then an 
extended QED Lagrangian includes the renormalisable terms 

n = 3 : fr075U M 7 / V (3.35) 
n = 4 : K pupa F^F p(7 + iipu p u u YD' y {c + d^)ip , (b,c,d : canst.) . 

We encountered the n = 4 terms before in Section 2: note that K^ vpa = 
U[^9u][ P Ua] has to obey rotation symmetry, which renders it equivalent 
to the term (I2.13p . and the last term above is equivalent to (12.111) (we 
proceed to covariant derivatives D u ). 

Refs. |189l |T90] explored the possibilities to add LIV operators of 
mass dimension n = 5, which are quadratic in the fields (of one kind) 
and not reducible to lower dimensions or total derivatives by applying 
the equations of motion. They assumed these (CPT odd) terms to be 
0(1/Mpi anc k) suppressedjff] In the pure gauge action, Ref. [189j identified 
only one (independent) term, 

( u p F pv {u ■ d)(u p F p ») , (F pv := \e^F aP ) . (3.36) 



^Planck 2 

The LIV parameter £ affects the photon dispersion relation as 

El = k 2 ±-^k 3 (k=\k\), (3.37) 

JWpianck 

where ± refers to left /right-handed circular polarisation^ 



58 0f course, extending this rule to lower dimensions n leads to another manifestation 
of the unsolved hierarchy problem, which we commented on before. 

59 Actually radiative corrections amplify the LIV term, leading to a fine-tuning prob- 
lem even for this dimension 5 term |191j . 

A perturbative quantisation for the corresponding scalar model — where the La- 
grangian involves second order time derivatives — was constructed in Ref. 192J. 
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The LIV modified fermion dispersion is similar, if we assume m <C 
E -C Mpi anc k (where m is the fermion mass). Then we can choose the 
spinor field to be an approximate eigenvector of 75 resp. of the chiral 
projectors, and we obtain [189J 

E 2 =p 2 + m 2 + V± -^- (p=\p\), (3.38) 

M Planck 

with r4 ermion = _^nti-fermion Ref discusses new q ED phenom- 

ena, which may occur due to this modification, including the transitions 

e~ — ► e~ + 7 vacuum Cerenkov radiation 

7 — > e + + e~ photon decay 

— > e + e + + e~ pair emission 

7 — > 27 , 37 ... photon splitting. (3.39) 

In Subsection 2.5 we already addressed vacuum Cerenkov radiation and 
the photon decay. In contrast, lepton pair emission and photon splitting 
cannot be arranged for by different MAVs, but by suitable SME param- 
eters [95l 1194] . Moreover we have discussed in Subsection 3.1 that the 
photon absorption 7 + 7 — > e + + e~ - which occurs in standard QED as 
well — may be altered by LIV terms, in particular in view of the energy 
threshold. 

A careful analysis has to keep track of angular momentum conserva- 
tion, which can be achieved e.g. in the photon decay above the energy 
threshold, since the momenta do not need to be collinear. Thus all the 
three LIV parameters, £, r] + and ^_ , are involved. This decay (as well 
as the pair emission) would occur rapidly above threshold, so that the 
observed high energy photons, in particular from the Crab Nebula and 
from blazars (with lower energy but longer time of flight, see Subsections 
3.1 and 3.2) set combined constraints on these three parameters [76J. 

If the GZK cutoff attenuates ultra high energy protons as predicted 
by standard field theory, this should also generate high energy photons, 

p + 7m -> A + -> p + 7T° -> p + 2 7uv . (3.40) 

Since it is not observed that UHECRs are accompanied by such UV pho- 
tons [23], Ref. |195] concludes that they are also subject to the expected 
attenuation, and derived a tiny bound for |£|, for the case rj± = 0. Ref. 
[196] generalised this study to obtain combined constraints for £ and rj±, 
considering in particular the processes 7 — > e + + e~ and 7 — > 37. 

Based on the Crab Nebula photons up to ~ 50 TeV one obtains 
extremely stringent constraints, if one relies on the standard assump- 
tion about the origin of its 7-radiation (synchrotron emission of high 
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energy electrons and positrons, and inverse Compton scattering of the 
synchrotron photons) [TBI f!97j. 



However, at this point we would like to be as economic as possible 
with our assumptions. Without any hypothesis about the mechanism 
causing the radiation, nor any assumption about the interactions in a 
LIV modified QED (which is needed for the transitions in the list (13.391) ) 
we can still address purely kinematic effects^ 

In this regard we have already discussed the energy dependence of the 
photon speed in Subsection 3.4. In addition the dispersion relation (13.371) 
also implies a helicity dependence for plane waves of the same momentum 
k. Over a distance I it leads to a relative delay Ath, 

t , ± = ^ 1± _^ _ At^MlL. (3.41) 

dk 2M Planck JWpianck 

It is a virtue of this formula that it is not affected by the uncertainty 
about a relative time-lag in the photon emission of different momenta k 
(unlike the studies reviewed in Subsection 3.4). However, one needs of 
course the (mild) assumption that there is no relative delay between the 
helicities at the source. 

As we mentioned before, Markarian 421 and 501 are the two known 
blazars which emitted photons that arrived with energies E > 10 TeV. 
The quasi-simultaneous arrival of both helicity states implies |193j 

|f | < 126 -> \M\ > 0.008 M Planck , (3.42) 

where M refers to the notation used in Subsection 3.4. This bound is in 
the same magnitude as those derived from the energy dependence. 

Birefringence : The helicity dependent energy (I3.37P yields yet an- 
other purely kinematic effect, which leads to a much stronger bound on 
the parameter f [76] @ In addition to the group velocity v± there are 
also distinct phase velocities. As a consequence, the polarisation plane 
of linearly polarised 7-radiation of momentum k is rotated after time t 
by the angle 

a(t) = i(E + - EJ) ^ . (3.43) 



60 Also in this respect there are attempts to go beyond the EFT approach and 
motivate the corresponding (parity breaking) modification of the Maxwell eqs. from 
quantum gravity [198j . 

61 Further studies address possible LIV birefringence effects in gravitational waves 
[199] and on the electron self-energy [63] . 
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Most useful for our purposes is the effect that a linearly polarised burst 
covering momenta k m { n . . . k max will be depolarised. The net polarisation 
is lost completely when the difference in the polarisation planes arrives 
at 

^ t / j 2 / 2 \ / \ 

QWx — a min = TT I "'max — "'mm) ~ 7T ' l"-^7 

iWpianck * 

We mentioned in Section 2 that this effect has been used to constrain 
the n = 3 Chern-Simons term ( 12. 9\\ [91], and it has also been applied to 
establish bounds on n = 4 LIV terms in the SME |200] . For the n = 5 
term that we are discussing here, this birefringence property was studied 
by Gleiser and Kozameh |131j . They referred to 7-rays originating from 
galaxies at a distance of 300 Mpc. For wave lengths ranging from 400 nm 
to 800 nm a linear polarisation was observed, within < 10° deviation of 
the polarisation plane. As a final highlight, let us quote that this yields 
a trans-Planckian bound of 

|e|<10- 4 - |M|>10 4 M Planck . (3.45) 

Of course this refers to a specific order in the EFT framework F 2 ! If 
we restrict the consideration to CPT even terms — returning to the 
concept outlined in Section 2 - - such operators of dimension n = 5 
are excluded. In the EFT scheme we would then naturally proceed to 
n = 6 operators causing LIV effects of 0(i? 2 /M| lanck ). For these terms 
helicity dependent photon velocities are still possible, but not required 
anymore [76] . Here the absence of vacuum Cerenkov radiation for protons 
of energies up to ~ 10 20 eV (see Subsection 2.5.2) becomes most powerful. 
As we anticipated in footnote[36l Ref. [103J discusses a link between SME 
parameters and hadronic effects by employing the parton model, though 
this connection is not ultimately clarified. 



4 Conclusions 

Cosmic rays involve particles at the highest energies in the Universe. 
Their flux is small at ultra high energies, which makes measurements 
difficult, but there is on-going progress in the techniques and extent of 
the corresponding observatories. In addition to the ultra high energy, 
we can also take advantage of their extremely long path to extract in- 
formation, which is inaccessible otherwise — that property is even more 
powerful in view of conclusions about scenarios for new physics. There- 
fore cosmic rays provide a unique opportunity to access phenomenology 

62 Returning to NC field theory, Ref. [201 j discusses how UV/IR mixing (see Sub- 
section 3.3) could affect birefringence. As for the laboratory, Ref. [202] proposed a 
high precision experiment on vacuum birefringence. 
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at tremendous energies. We have seen in the last Subsections that this 
can reach out to energies only two or three orders of magnitude below the 
Planck scale, and for specific parameters one even attains trans-Planckian 
magnitudes. Hence predictions in this regime have now arrived at falsi- 
fiability, which is a basic criterion for a scientific theory [203J. Although 
this only leaves a small window for possible LIV parameters, there is a 
wide-spread consensus that this scenario deserves attention — over the 
past 10 years a preprint on this subject appeared on average about once 
a week, and at least 52 are top-cited. 

We have discussed the challenges to verify the GZK cutoff for UHE- 
CRs, and the multi-TeV 7-puzzle. The question is: why is the Universe 
- with its Cosmic Background Radiation (CBR) — so transparent for 
protons, or heavier nuclei, with energies E p > 6 • 10 19 eV, and for photons 
with Ery > 10 TeV ? In both cases it is still controversial if the puzzle 
really persists. 

The photopion production of UHECR protons with CMB photons is 
actually a process of harmlessly low energy (around 200 MeV) in the 
rest frame of the proton. It can be studied in detail in our laboratories. 
However, the application of the results for the cross-section etc. to UHE- 
CRs requires boosts with Lorentz factors of the order 7 ~ O(10 n ). The 
applicability of such extreme Lorentz transformations cannot be tested 
in the laboratories (which are limited to 7 < O(10 5 )). Hence the viability 
of LI under extreme boosts is a central issue. If a puzzle persists at last, 
LIV provides a possible solution. In Section 2 we discussed a theoretical 
basis of LIVs (SME and in particular MAVs), and its impact on cosmic 
rays physics — a tiny LIV could remove the GZK cutoff. 

As an alternative theoretical framework for LIV, Section 3 addressed 
field theory in a NC space. This approach is fundamentally different 
from the effective action concepts that we considered in Section 2. We 
focused in particular on the dispersion relation of the photon in a NC 
world. A possible distortion of the photon dispersion can be probed 
best by the observation of GRBs or blazar flares. We also discussed 
the Effective Field Theory approach to LIV terms for cosmic photons 
and stringent limits on leading terms, based on momentum or helicity 
dependent photon velocities. 

Up to now, no experimental evidence for LIV has been found any- 
where, despite more and more precise tests. As examples we reviewed 
attempts with cosmic neutrinos and photons (other high precision tests 
are performed in atomic physics [721 [TTJ EH]). However, a large number 
of LIV parameters are conceivable. For some of them the tested precision 
may still be insufficient to rule out interesting scenarios for new physics, 
which could shed completely new light on the open issues mentioned 
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above. If the GZK cutoff is indeed confirmed — as the data by the col- 
laborations HiRes and Pierre Auger suggest — then this could be viewed 
as indirect evidence in favour of the validity of LI even at 7 ~ O(10 n ). 

New projects to detect UHECRs are on the way. As examples we 
mention Telescope Array (TA) |204] (another hybrid approach with an 
air shower array on ground and fluorescence telescopes), the Extreme 
Universe Space Observatory (EUSO) [205] and the Orbiting Wide-angle 
Light-collectors (OWL) [206j (both are going to search for air fluorescence 
light from satellites). The general prospects for observatories in space are 
discussed in Ref. |2U7j . 

High energy j-rays are detected by satellite-based telescopes in the 
AGILE Observatory [208] (Astro-rivelatore Gamma a Immagini Leg- 
gero), and in the GLAST Observatory [185] (systematic observation of 
cosmic photons over a broad energy). On ground there are suitable 
Cerenkov detector arrays, such as the High Energy Stereoscopic System 
H.E.S.S. in Namibia [IE] and MAGIC (we have already referred to 
GLAST and MAGIC in Section 3). Ref. |21U] gives a recent overview. Re- 
garding GRBs, the Chinese-French Space-based multi-band astronomical 
Variable Objects Monitor mission (SVOM) [211] plans to detect about 
80 bursts a year. 

As for the upcoming search for high energy cosmic neutrinos (out- 
lined e.g. in Ref. [128]) we mention the IceCube Neutrino Observatory 
|212] . which detects Cerenkov light in the antarctic ice, as well as ANITA 
(Antarctic Impulse Transient Antenna, balloon-borne), ANTARES (As- 
tronomy with a Neutrino Telescope and Abyss environmental RESearch, 
at sites deep in the Mediterranean Sea) and NESTOR (Neutrino Ex- 
tended Submarine Telescope with Oceanographic Research) |213j . 

Let us also mention a completely new technique, which is currently 
under investigation |214] : the idea is to measure the synchrotron radia- 
tion of air shower leptons in the magnetic field of the Earth, in particular 
radio flashes, with an array of 30 dipole antennae. 

Finally we should emphasise once more that the Pierre Auger Ob- 
servatory is still at an early stage; it is operating successfully since the 
beginning of the year 2004, but the final Cerenkov detectors have been 
installed in 2008. Moreover that Collaboration plans another observatory 
in the northern hemisphere [12]. 

This is a very active field of research with exciting open questions. 
We may expect enlightening new data in the near future. They could 
lead to new insight in outstanding issues like LIV — or to new puzzles 
and perhaps to evidence for new physics. 
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A New development since Nov. 2007: 
The AGN Hypothesis 

In November 2007 the Pierre Auger Collaboration published a spectac- 
ular report |215j . which is relevant for the issues that we discussed in 
Section 1. They analysed the arrival directions of the UHECRs that 
they had detected so far. We denote their observation as the AGN Hy- 
pothesis: 

The UHECRs directions are clustered and correlated with the locations of 
nearby Active Galactic Nuclei (AGN). 

This is clearly in contrast to the traditional assumption of an isotropic 
distribution (see Section 1.1)1^1 The possibility that the AGN — or other 
objects located next to them — could be the UHECR sources appears 
consistent because they refer to AGN nearby: then the UHECR may 
be little deflected by magnetic fields, and attenuated only mildly by the 
CMB. In this scenario one could assume also AGN far away to emit 
UHECRs, but they would arrive with energies below the GZK cutoff 
at Earth. Therefore the AGN Hypothesis agrees with the validity of 
the GZK cutoff; it disfavours the possibility that UHECRs indicate new 
physics. In particular Ref. [218] deduces an upper bound for the MAV 
difference cv — c p < 5 • 10~ 23 . On the other hand, it could be the starting 
point of a new epoch of cosmic ray astronomy. 

The Pierre Auger Collaboration started from the UHECR data col- 
lected in the period from January 2004 to May 2006. Their analysis 
involved 3 parameters 

E m[n : threshold energy to be counted as an UHECR 

■0 : angular interval around an UHECR direction 

i? max : maximal distance to a "nearby AGN" . 
Thanks to the analysis of hybrid events, which were detected by the 
extended air shower and by the fluorescence telescopes, the energy has 
a relatively small uncertainty around 22% (cf. Subsection 1.4), and the 
angular accuracy is about 1° [11]. E min and -R max were of course assumed 
in the order of the GZK cutoff and the predicted super-GZK path length 
(eqs. (II. 2p and (j!.15p ). The distance to the AGN is known from the 
redshift of their 7-radiation. 

63 For the characteristic of AGN we refer to footnote [71 and for an earlier suggestion 
that they could be the UHECR source to Ref. |216| . The question if they could 
be endowed with a sufficiently powerful acceleration mechanism has recently been 
reconsidered in Refs. [217] . 



62 



These 3 parameters were tuned to the values 



(£ min , V, #ma X ) = (5.6 ■ 10 19 eV, 3.1°, 75 Mpc) , (A.l) 

which gave maximal support for the AGN Hypothesis: it captures 12 
out of 15 UHECRs, whereas only 3.2 would be expected for hypothetical 
isotropic rays (although in that case one could question if one should 
keep these three parameters fixed). 

As a clean test this Collaboration then analysed — still with fixed 
parameters flA.lj) — the UHECRs that they detected later, in the period 



May 2006 to August 2007. In that case it captures 8 out of 13 UHECRs. 
So the ratio went down, but the number 8 is still large compared to 2.8, 
which would be the expectation for isotropic rays. 

The ensemble of UHECR directions observed by the Pierre Auger 
Collaboration in the whole period January 2004 to August 2007 is shown 
in the Figure H5J it covers their 27 top UHECRs (E > 5.7 • 10 19 eV). 
The black circles mark the range within 3.1° of their directions of arrival, 
and the asterisks correspond to the 472 known nearby AGN (up to the 
distance -R max in eq. flA.lj) . which corresponds to redshift z = 0.018), 
according to the Veron-Cetty & Veron Catalogue |219j . The alternative 
use of the Swift BAT Catalogue leads qualitatively to the same results 
In further investigations Ref. [221] (Ref. |222j ) reports a correlation 



between these 27 UHECR directions and the locations of spiral galaxies 
(extended nearby radiogalaxies) by using the HIPASS Catalogue (several 
catalogues). 



A.l Comments on the AGN Hypothesis 

The AGN Hypothesis consists actually of two parts, and we consider it 
worthwhile discussing them separately: 

(a) Clustering of the UHECR directions. 

(b) Correlation of the clustered directions with locations of nearby AGN. 
Part (b) depends on (a), but the validity of (a) only could be an option. 

Ref. |215j reports 99 % C.L. for part (a) and here only two tuned 
parameters enter, E m [ n and ip. Still this new UHECR statistics is in 
the same magnitude as the total previous world data (and the same 
holds for the total exposure so far, i.e. [detection time] x [area]; a com- 
parative Table is given in Ref. [52]). Hence one may wonder why this 
clustering did not become obvious earlier. Of course, previous data were 
also analysed with this respect, and some indication for clustering was 
found in the world data available in 1995 (from Volcano Ranch, Hav- 
erah Park, Yakutsk and AGASA at an early stage) [223]. This inspired 
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Figure 15: A projection of the celestial sphere (plot adapted from Ref. 
]21h^ ): the circles mark the arrival directions of the 27 UHECRs detected 
by the Pierre Auger Collaboration, and the asterisks the locations of 472 
known nearby AGN. In particular for Centaurus A a white symbol is 
used. The dashed line follows the super- galactic plane. The darkness of 
the areas corresponds to the exposure (in the white area it vanishes). 

the AGASA Collaboration to perform a careful analysis of the directions 
where they found UHECRs [224J. They could see a slight clustering 
near the super-galactic planer 4 ! but it was weaker than in the previous 
world data; otherwise they found consistence with isotropy. Not even 
their slight anisotropy signal was confirmed by HiRes, which reported 
full isotropy |225j . A new analysis referring to an updated definition of 
the super-galactic plane confirms a clustering of the 27 top UHECR di- 
rections observed by the Pierre Auger Collaboration, but for a similar 
number of previous world data events (in the same energy magnitude) 
clustering is still not evident [226J. 

With respect to part (b) of the AGN Hypothesis the statements about 
the C.L. are less striking in Ref. [215j . The analysis in Ref. [227] is 
more complete, for instance regarding the variation of the parameters 
(lA.lj) (see also Ref. |228j ). First support for this Hypothesis beyond the 
Pierre Auger Observatory was provided by a new analysis of the Yakutsk 
data |229j . but not by the HiRes data [230 : in that case, only 2 out 
of 13 events fulfil the correlation criterion for the parameter set ( lA.lj) . 

64 The super-galactic plane is a sheet-like structure that contains the local super- 
clusters. It corresponds to the dash line in Figure [T51 
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Ref. |231j comments that the effect could reflect the property that the 
known AGN, as well as the detected nearby UHECR sources, are both 
preferentially located near the super-galactic plane. 

Critics was also expressed in Ref. [232] , which starts from the state- 



ment that the arriving flux of a source at distance R should be oc 1/R 2 . 
Hence according to the AGN Hypothesis the nearest AGN should con- 
tribute most. Therefore Ref. [232] expects Centaurus A and Virgo - 
which contain the two nearest AGN — to contribute each about 6 out 
of the 27 UHECR in the sky-map of Figure HH if t ne AGN Hypothesis 
holds. The direction of the galaxy Centaurus A does contribute, and 
(with some tolerance about the angular window) it is compatible with 
this expectation^ but nothing is found in the direction of the Virgo 
cluster of galaxies (in Figure [15] it is located in the upper part near 
the dashed line, which follows the super-galactic plane). That deficit is 
the main point in the analysis of Ref. [232] . which claims that the data 
disfavour the AGN Hypothesis with 99 % C.L. One might remark that 
this reasoning seems to refer to a scenario where all AGN are UHECR 
sources with continuous luminosity. If AGN are episodic sources, then 
this argument might be circumvented!^! 

For sure the matter distribution in the GZK-sphere is inhomogeneous. 
Therefore the AGN hypothesis could indeed lead to a new branch of 
astronomy based on cosmic rays, as we mentioned in the beginning of 
this appendix. A possible correlation of the nearby matter distribution 
with the magnetic field strength is another controversial issue [236] . The 
way it could affect the UHECR propagation is discussed — for various 
scenarios of the magnetic field structure — in Ref. [12] . 

Finally we add another interpretation put forward in Ref. [237] , which 
concludes from the UHECR path limitation — according to the analysis 
in Ref. [215J - that the primary particles at source could be mostly 
nuclei heavier than protons. In fact the distance R mSiX to the nearby 
AGN is a bit shorter than required for the maximal super-GZK path 
length of protons (eq. (11.151) ). hence better suited for other nuclei. A new 
calculation of the effective energy loss of UHECRs in Ref. [231] supports 
the conjecture that these rays are mainly emitted as light nuclei, up to 
an atomic mass number A « 24. This interpretation caused a surge of 
interest in the CBR attenuation of nuclei, which has also been revisited in 
Refs. [2381 HE] @ Ref. [239] added that this could affect the magnitude to 



65 See Ref. |233| 1234] for discussions of the high energy radiation mechanism of 
Centaurus A. 

66 Ref. [235] compares the scenarios of continuous vs. episodic UHECR sources, and 
arrives at different predictions for the energy dependence of the 7 and v flux. 

67 One expects UHECRs to be accompanied by a neutrino flux generated in CMB 
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which LIV parameters are ruled out if the AGN Hypothesis holds, since 
possible LIV effects in the dispersion of N bound particles are amplified 
oc iV 2 . On the other hand, for a given energy the Lorentz factor between 
the FRW laboratory frame and the ray rest frame shrinks oc 1/A. 

After a long journey, Ref. [240] expects either protons or Fe nuclei (or 
perhaps sub-Fe nuclei) to survive as UHECR primary particle. In this 
regard, a statistical analysis of the penetration depth into the atmosphere 
(cf. eq. (I1.17P ) suggests a mixed composition |241j . However, also that 
observation appears to be consistent only in part with the HiRes data 
[242J (a comparison is shown in Figure 4 of Ref. [243J ) . 

In any case, more statistics will be needed to identify the nature of 
the UHECR primary particles directly. Above 10 19 eV not much this is 
known in this respect. We may hope for new insight based on upcoming 
data, which the Pierre Auger Collaboration intends to release in 2009. 

B Table of short-hand notations 

For convenience we add an alphabetic list of the short-hand notations 
used in the text, with some hints where the corresponding terms are 
introduced and used: 

AGN Active Galactic Nuclei, see footnote [TJ 
They play a central role in Appendix A. 

C.L. Confidence Level (for statistical results). 

CPT Charge conjugation combined with parity and time reflexion. 
Subsection 2.1 comments on its relation to LIV. 

CMB Cosmic Microwave Background, see Subsection 1.2. 

DSL Double Scaling Limit, see Subsection 3.3.2. 

DSR Doubly Special Relativity, see Subsection 3.3. 

EFT Effective Field Theory, see Subsection 3.5. 

FRW Friedmann-Robertson- Walker 

(or Friedmann-Lemaitre- Robertson- Walker) 

A line element in the FRW metrics is given by 

ds 2 = dt 2 — a(t) 2 (dr 2 + f 2 dVt) (dVl: spherical element) 

a(t): scale factor of the Universe, f : covariant distance 

f = r (no curvature), f = R sin {r/R) (f = Rsinh(r/ R)) 

for positive (negative) curvature with radius R. 

GRB Gamma Ray Burst, see Subsection 3.2. 

GRT General Relativity Theory 

GZK Greisen-Zatsepin-Kuz'min 

interactions (see e.g. Ref. [128] ). which also depends on the atomic mass A [234\. 



66 



The GZK cutoff for the energy of cosmic rays, 
Eqzk — 6 • 10 19 eV, is discussed in Subsection 1.3. 

LI Lorentz Invariance (or Lorentz Invariant) 

LIV Lorentz Invariance Violation, addressed in Sections 2 and 3. 

MAV Maximal Attainable Velocity, introduced in Subsection 2.4. 

NC Non-Commutative, see Subsection 3.3. 

SSB Spontaneous Symmetry Breaking 

SME Standard Model Extension, see Subsection 2.3. 

SRT Special Relativity Theory 

SUSY Supersymmetry 

UHECR Ultra High Energy Cosmic Ray, a cosmic ray with energy 
near (or above) the magnitude of E GZ k 

For completeness we add that the literature often expresses the huge 
energies involved in cosmic rays in units of PeV, EeV and ZeV, where 
the initial capital letters mean Peta : 10 15 , Exa : 10 18 , Zetta : 10 21 . 
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